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INTRODUCTION 

The  problem  of  characterizing  shock  wave  properties  in  solid  materials 
necessitates  a  knowledge  of  the  solid  equation  of  states  at  high  pressures. 

The  pressures  attainable  in  strong  shock  waves  can  exceed  considerably 
the  pressures  attainable  by  experimental,  static  techniques.  Static  pressures 
which  have  been  attained  experimentally  are  about  in  the  hundred  kilobar  range, 
whereas  dynamically,  pressures  of  about  five  megabars  have  been  attained. 
Thus,  it  has  been  necessary  to  use  dynamic  techniques  to  determine  equations 
of  state  for  use  in  shock  wave  calculations. 

One  common  technique  to  obtain  an  equation  of  state  at  high  pressure  is 
to  calibrate  an  assumed  equation  of  state,  i.e.,  the  Mie-Gruneiscn  equation  of 
state,  in  terms  of  a  shock  Hugoniot.  The  Mie-Griineisen  equation  of  state  is 

P  -  P„  *  $  (E  -  E„) 

where  Y  is  the  Griineisen  parameter.  For  the  state  of  a  material  after  shock 
passage,  one  can  write 

PH  '  P0  =  V  *EH  '  E0> 


Combining  the  above  two  equations  and  assuming  that  the  Gruneisen 
parameter  is  a  function  of  volume  only  gives 

p  ■  pn  '  v  <E  -  V 

The  above  equation  represents  an  equation  of  state  for  small  excursions 
from  the  calibrating  Hugoniot.  I  to  validity  depends  on  the  assumptions  that 
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the  compression  at  is  the  same  as  would  be  produced  by  a  hydrostatic 
pressure  of  the  same  magnitude;  the  Hugoniot  state  after  shock  passage  is 
an  equilibrium  state;  and  that  the  Gruneisen  parameter  does  not  vary  with 
temperature  but  is  a  function  of  volume  only. 


Recently,  shock  work  with  foams  and  measurements  of  the  sound 
velocity  behind  shock  waves  have  raised  the  question  as  to  whether  Y  is  indeed 
independent  of  temperature. 

Also  a  problem  in  calculating  the  shock  properties  of  foamed  and  dis¬ 
tended  materials  is  how  to  formulate  equations  of  state  in  the  lower  pressure 
regions.  At  lower  pressures,  the  assumption  of  these  materials  attaining 
states  of  "true"  thermodynamic  equilibrium  may  no  longer  be  tenable.  (The 
term  true  is  put  in  quotes  because,  as  is  discussed  in  Section  6B,  the  defini¬ 
tions  of  various  equilibria  are  rather  confused.) 

In  trying  to  evaluate  the  possible  temperature  dependency  of  Y,  it  was 
soon  found  that  the  topic  of  the  Gr&neisen  parameter  is  treated  in  the  litera¬ 
ture  in  a  limited,  fragmentary  way.  Nowhere  floes  there  exist  a  systematic, 
overall  treatment.  Also,  there  appears  to  be  a  confusion  about  possible 
definitions  of  the  Gruneisen  parameter.  Thus,'  the  task  objective,  as  pre¬ 
sented  in  Experimental  Memo  FR961,  set  for  the  present  study  was;’ 

“Contribute  to  better  understanding  of  the  Gruneisen  parameter}-* 
its  theoretical  basis,  experimental  determination,  and  usage  in  present 
hydrodynamic  response  code.** 

The  tasks  set  for  this  study  were; 

Task  1  -  Compile  systematic,  overall  treatment  of  theoretical 
t  aspects  of  the  Gruneisen  parameter  (starting  with  original 
Grflneisen  work,  circa  1912).  Include  complete  bibliography 
and  pertinent  excerpt  from  source  documents. 

Task  2  -  Review  aH  experimental  techniques  (past  and  present)* 
used  to  derive  experimental  values  of  Gruneisen  parameter. 

Compare  data  for  selected  materials  (choose  materials  relatively 
insensitive  to  fabrication  variables  and  anisotropic  effects). 
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Suggest  most  useful  approaches  for  typical  spectrum  of  reentry 
vehicles  materials,  i.e.,  organic  composites,  epoxy  adhesives 
and  light  metals.  Include  complete  bibliography  and  pertinent 
excerpts  from  source  documents.  Describe  experimental  setups 
in  detail. 

Task  3  -  Calculate  Grfineisen  parameter  for  aluminum,  beryllium, 
tungsten,  tantalum,  tungsten  carbide  and  silicon  carbide  as  function 
of  pressure  to  30  kilobars  and  temperature  up  to  melt.  Do  same  for 
nylon  phenolic  and  refrasil  phenolic,  to  50  kilobars  pressure  and 
temperatures  of  500°F. 


During  the  last  three  months  of  the  study  work  was  concentrated  on  the 
following  task,  a's  presented  in  the  Task  Memo  TP 7087: 

"Continue  the  Gruneisen  parameter  study  in  the  area  of  porous, 
foamed  and  distended  materials.  Review  all  techniques  used  for  deriving 
the  Gruneisen  parameter  for  such  materials  and  present  any  pertinent  data. 
Make  recommendations  for  the  incorporation  of  this  parameter  into  our 
present  computer  code." 


The  preceding  tasks  are  covered  in  this  report  in  the  following  manner. 
Fundamental  thermodynamic  and  physical  definitions  are  covered  in  Section  1. 
The  theoretical  bases  of  equations  of  state  and  the  Mie -Gruneisen  Equation 
of  State  are  discussed  in  Sections  2  and  3.  The  experimental  techniques  of 
determining  the  various  Gruneisen  parameters  for  solid  materials  are 
reviewed  in  Section  4,  and  actual  calculations  of  the  various  Gruneisen 
parameters  for  the  materials  of  the  study  are  presented  in  Section  5.  The 
thermodynamics  of  foamed  and  distended  materials  are  considered  in 
Section  6. 
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FOREWORD 

The  study,  performed  at  Lockheed's  Huntsville  Research 
it  Engineering  Center  during  the  period  28  September  1966  to 
27  October  1967,  is  in  accordance  with  the  task  objectives  and 
task  statements  presented  in  the  Lor kheed/Sunnyvale  Experi¬ 
mental  Memos  PR961  and  TP70R7.  Technical  direction  was 
furnished  by  Mr.  D.  Fenton  (Dept.  8!  -21)  and  Mr.  E.  L.  Esch 
(Dept.  31-73).  Their  support  and  encouragement  are  gratefully 
acknowledged. 

Mr.  Ivan  Landis  contributed  significantly  to  the  study  by 
performing  the  mathematical  analysis  reported  in  Section  SC  and 
conducting  the  study  phase  reported  in  Sections  5D  and  SE. 

Mr.  Ernest  Raper  performed  the  calculations  presented  in 
Section  5E. 
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SUMMAR  Y 

The  theoretical  basis  of  the  Mi  Griineisen  equation  of  state  was  con¬ 
sidered  and  systematized.  The  most  important  experimental  techniques  for 
determining  the  several  Griineisen  parameters  were  reviewed.  Several 
Griineisen  parameters  for  some  of  the  materials  of  the  study  were  calculated 
hy  methods  presented  in  the  literature  and  by  an  original  method  developed 
in  the  course  of  the  study.  Also,  the  prohlems  involved  in  characterizing  the 
states  of  foamed  materials  before  and  after  irreversible  changes  were 
reviewed. 

The  most  important  facts  drawn  from  the  study  of  the  theoretical  basis 
of  the  Mie-Griineisen  equation  of  state  are  (a)  the  Mie-Grunei Ben  equation  is 
necessarily  valid  if  y^,  not  y^,  in  a  certain  temperature  range  is  independent 
of  temperature  and  (b)  the  Mie-Griineisen  equation  is  necessarily  valid  if  Y^ 
is  temperature  independent  throughout  the  temperature  range  from  0°K 
upwards. 

Although  our  calculations  of  y,  at  low  pressures  for  the  metals  of  the 
study  do  not  show  much  temperature  variation  from  room  temperature  up, 
values  of  Y^  it  very  low  temperatures  found  in  the  literature  for  aluminum  do 
vary  considerably  as  apparently  do  the  y^  values,  which  we  have  calculated  by 
the  method  of  Beecroft  and  Swenson  and  our  own  method.  Thus,  it  would 
appear  that  the  use  of  the  Mie-Griineisen  equation  of  state  for  metals  at  low 
temperatures  is  not  valid. 

At  moderate  pressures  our  calculations  of  Y ^  for  sodium  and  aluminum 
indicate  that  the  variation  of  Y ^  with  temperature  is  about  the  same  as 
observed  at  low  pressures  with  perhaps  very  unusual  behavior  in  the  very 
high  pressure  regions.  However,  time  limitations  did  not  permit  a  check  of 
the  sensitivity  of  the  calculation  to  the  approximations  made,  so  the  results 
cannot  be  considered  conclusive. 
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The  conclusions  derived  from  a  review  of  (he  thermodynamics  of 
foamed  and  distended  materials  are  that  considerable  confusion  exists  ahout 
basic  definitions  of  various  types  of  equilibrium  and,  therefore,  about  calcu¬ 
lation  of  thermodynamic  state  functions  for  materials  which  have  undergone 
changes  involving  hysteresis  of  one  sort  or  another.  However,  a  classical 
thermodynamic  approach  to  "totally''  irreversible  processes  promises  to 
yield  much  information  once  the  discussions  of  various  authors  are  clarified 
and  correlated.  The  preceding  state  of  affairs  is  really  pertinent  to  the  cal¬ 
culations  of  a  Gruneiscn  parameter  in  the  low  pressure  region.  The  high 
pressure  region  where  the  foam  is  assumed  completely  crushed  and  annealed 
can  be  handled  straightforwardly. 

Time  limitations  did  not  permit  the  full  exploration  of  the  effects  of 
material  anisotropy,  material  rigidity,  melting,  or  material  porosity  and 
inhomogeneity.  Because  the  problems  associated  with  the  last  two  factors 
have  not  yet  been  resolved.  Griineisen  parameters  for  the  plasticB  of  the 
study  were  not  calculated. 
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see  p.  2-7 
see  p.  3-2 
see  p.  4-8 
see  p.  6-4 
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INTRODUCTION 

The  problem  of  characterizing  shock  wave  properties  in  solid  materials 
necessitates  a  knowledge  of  the  solid  equation  of  states  at  high  pressures. 

The  pressures  attainable  in  strong  shock  waves  can  exceed  considerably 
the  pressures  attainable  by  experimental,  static  techniques.  Static  pressures 
which  have  been  attained  experimentally  are  about  in  the  hundred  kilobar  range, 
whereas  dynamically,  pressures  of  about  five  megabarB  have  been  attained. 
Thus,  it  has  been  necessary  to  use  dynamic  techniques  to  determine  equations 
of  state  for  use  in  shock  wave  calculations. 


One  common  technique  to  obtain  an  equation  of  state  at  high  pressure  is 
to  calibrate  an  assumed  equation  of  state,  i.e.,  the  Mie-Griineisen  equation  of 
state,  in  terms  of  a  shock  Hugoniot.  The  Mie-Griineisen  equation  of  state  is 

P  -  P0  =  £  (F.  -  E#) 


where  Y  is  the  Griineisen  parameter.  For  the  state  of  a  material  after  shock 
passage,  one  can  write 


P 


H 


Y 

V 


<eh  -  Eo> 


Combining  the  above  two  equations  and  assuming  that  the  Griineisen 
parameter  is  a  function  of  volume  only  gives 

P  -PH  =  $<e.-eh) 

The  above  equation  represents  an  equation  of  state  for  small  excursions 
from  the  calibrating  Hugoniot.  Its  validity  depends  on  the  assumptions  that 
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the  compression  at  is  the  same  as  would  be  produced  by  a  hydrostatic 
pressure  of  the  same  magnitude;  the  Hugoniot  state  after  shock  passage  is 
an  equilibrium  state;  and  that  the  Gruneisen  parameter  does  not  vary  with 
temperature  but  is  a  function  of  volume  only. 


Recently,  shock  work  with  foams  and  measurements  of  the  sound 
velocity  behind  shock  waves  have  raised  the  question  as  to  whether  Y  is  indeed 
independent  of  temperature. 

Also  a  problem  in  calculating  the  shock  properties  of  foamed  and  dis¬ 
tended  materials  is  how  to  formulate  equations  of  state  in  the  lower  pressure 
regions.  At  lower  pressures,  the  assumption  of  these  materials  attaining 
states  of  "true"  thermodynamic  equilibrium  may  no  longer  be  tenable.  (The 
term  true  is  put  in  quotes  because,  as  is  discussed  in  Section  6B,  the  defini¬ 
tions  of  various  equilibria  are  rather  confused.) 

In  trying  to  evaluate  the  possible  temperature  dependency  of  Y,  it  was 
soon  found  that  the  topic  of  the  Gruneisen  parameter  is  treated  in  the  litera¬ 
ture  in  a  limited,  fragmentary  way.  Nowhere  does  there  exist  a  systematic, 
overall  treatment.  Also,  there  appears  to  be  a  confusion  about  possible 
definitions  of  the  Gruneisen  parameter.  Thus,  the  task  objective,  as  pre¬ 
sented  in  Experimental  Memo  PR961,  set  for  the  present  study  was; 

"Contribute  to  better  understanding  of  the  Gruneisen  parameter;  • 
its  theoretical  basis,  experimental  determination,  and  usage  in  present 
hydrodynamic  response  code." 

The  tasks  set  for  this  study  were: 

Task  1  ••  Compile  systematic,  overall  treatment  of  theoretical 
,  aspects  of  the  Gruneisen  parameter  (starting  with  original 
Gruneisen  work,  circa  1912).  Include  complete  bibliography 
and  pertinent  excerpt  from  source  documents. 

Task  2  Review  all  experimental  techniques  (past  and  present) 
used  to  derive  experimental  values  of  Gruneisen  parameter. 

Compare  data  for  selected  materials  (choose  materials  relatively 
insensitive  to  fabrication  variables  and  anisotropic  effects). 
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Suggest  most  useful  approaches  for  typical  spectrum  of  reentry 
vehicles  materials,  i.e.,  organic  composites,  epoxy  adhesives 
and  light  metals.  Include  complete  bibliography  and  pertinent 
excerpts  from  source  documents.  Describe  experimental  setups 
in  detail. 

Task  3  -  Calculate  Griineisen  parameter  for  aluminum,  beryllium, 
Tungsten,  tantalum,  tungsten  carbide  and  silicon  carbide  as  function 
of  pressure  to  30  kilohars  and  temperature  up  to  melt.  Do  same  for 
nylon  phenolic  and  refrasil  phenolic  to  SO  kilobars  pressure  and 
temperatures  of  500°F. 


During  the  last  three  months  of  the  study  work  was  concentrated  on  the 
following  task,  as  presented  in  the  Task  Memo  TP  7087: 


"Continue  the  Griineisen  parameter  study  in  the  area  of  porous, 
foamed  and  distended  materials.  Review  all  techniques  used  for  deriving 
the  Griineisen  parameter  for  such  materials  and  present  any  pertinent  data. 
Make  recommendations  for  the  incorporation  of  this  parameter  into  our 
present  computer  code." 


The  preceding  tasks  are  covered  in  this  report  in  the  following  manner. 
Fundamental  thermodynamic  and  physical  definitions  are  covered  in  Section  1. 
The  theoretical  bases  of  equations  of  state  and  the  Mie-Grvineisen  Equation 
of  State  are  discussed  in  Sections  2  and  3.  The  experimental  techniques  of 
determining  the  various  Griineisen  parameters  for  solid  materials  are 
reviewed  in  Section  4,  and  actual  calculations  of  the  various  Griineisen 
parameters  for  the  materials  of  the  study  are  presented  in  Section  5.  The 
thermodynamics  of  foamed  and  distended  materials  are  considered  in 
Section  6. 
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Section  1 

ELEMENTARY  PHYSICS  AND  THERMODYNAMICS  OF  SOLIDS 

The  total  internal  energy  of  any  isolated  system  is  made  up  of  a  kinetic 
part  and  a  potential  part. 


E  =  K  +  U  (1) 

The  potential  energy  of  a  solid  has  its  origin  in  the  attractive  and  re¬ 
pulsive  forces  between  the  atoms,  molecules,  or  ions.  In  an  ideal  gas  it  is 
assumed  that  the  forces  between  the  molecules  are  so  small  that  the  average 
potential  energy  is  negligible  compared  to  the  average  kinetic  energy.  In  a 
solid,  on  the  other  hand,  the  potential  energy  plays  a  very  important  role. 

The  internal  energy  of  a  solid  is  often  represented  as  (Fumi  and  Tosi,  1962) 

E(V,  T)  =  Eg( V)  +  F.V(V,T)  (2) 

where  Eg(V)  represents  purely  volume -dependent  static  energy  of  the  solid, 
or  lattice  energy,  and  EV(V,  T)  represents  the  vibrational  energy.  The  term 
Eg  can  be  further  defined  as  the  potential  energy  if  the  amplitude  of  the  ther¬ 
mal  vibrations  were  all  zero.  Fletcher  (1997)  points  out  that  in  general  Eg 
is  not  a  function  of  volume  only.  However,  it  is  a  fair  approximation  where 
electronic  effects  are  negligible.  Thus,  we  will  not  concern  ourselves 
further  with  the  possible  temperature  dependence  of  Eg. 

Another  equivalent  way  of  representing  the  energy  is 

E(V,  T)  =  Eq(V)  +  ET(V,  T)  (3) 

where  Eq(V)  is  the  energy  of  the  solid  at  absolute  zero,  and  E^,(V,  T)  repre¬ 
sents  the  thermal  energy. 
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At  zero  degrees  absolute,  the  term  E^.(V,  T)  becomes  zero,  but  Ey(V,  T) 
does  not.  The  term  Ey(V,T),  at  absolute  zero,  is  designated  as  the  zero- 
point  vibrational  energy.  Thus,  the  relation  between  Eg  and  Eg  can  be  written 

Eq(V)  =  ES(V)  +  E°(V)  (4) 

where  Ey  is  a  function  of  volume  only  and  represents  the  zero-point  vibrational 
energy.  In  statistical  mechanics  Ey  is  given  as  (Slater,  1939,  p.  216) 

3N 

EV  l/2(hi/.>  (5) 

j=l 


where  V.  is  the  frequency  of  the  j  mode  of  oscillation. 
J 


In  materials  where  the  interatomic  forces  and  the  atomic  masses  are 
great,  the  zero-point  energy  is  small  and  can  usually  be  neglected.  However, 
in  a  substance  like  helium,  the  zero-point  energy  iB  very  significant.  (Zemansky, 
1957,  p.  378) 


The  general  variation  of  Eg  as  a  function  of  volume  is  shown  in  Figure  1-1. 


Figure  1-1  -  Energy  and  Pressure  as  a  Function  of  Volume  at  Zero  Degrees 
Absolute  (after  Zel'dovich,  Ya.  B,  and  Yu.  P.  Rayzer,  1965*  p.  843) 
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As  Figure  1-1  shows,  Eg  is  a  minimum  atV^^,  the  volume  at  zero 
degrees  absolute  and  zero  pressure.  (NOTE:  The  energy  at  this  point  Egg 
is  arbitrarily  assigned  the  value  of  zero).  The  pressure,  Pg,  at  absolute 
zero  is  sometimes  called  the  elastic  pressure  and  is  given  by: 


(6) 


The  term  Eg  for  volumes  not  too  far  removed  from  Vgg,  can  be  ex¬ 
panded  to  (Mott  and  Jones,  1958,  p.  16) 


The  second  term  on  the  right  vanishes  because  d  Eg/d  V  at  Vgg  is  zero. 

The  term  Eg  can  be  similarly  expanded,  but  the  volume  at  which 
dEg/dV  is  zero  is  slightly  different  from  Vgg  (Huang,  1951,  p.  205).  In 
this  case  the  volume  should  be  tne  static  equilibrium  volume.  This  dis¬ 
tinction  appears  to  be  generally  ignored  or  confused  in  the  literature. 


Using  the  first  or  vibrational  formulation  the  energies  can  be  further 
broken  down  to 


E(V ,  T)  =  Ky(V,T)  +  Uy(V,T)  +  Eg(V)  (8) 

The  terms  Ky  and  Uy  represent  the  kinetic  and  potential  parts  of  the  vib¬ 
rational  energy.  In  statistical  mechanics  equation  8  is  written  (Slater,  19-?, 
p.  21 7) 

3N 

E  =  E  (V)  +]P  n  h 


(9) 
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where  n.  represents  the  quantum  number  of  the  j  vibration  mode.  The  term 
J 

3N 

n.  h  v.  includes  Ky  and  Uy,  that  is,  the  vibrational  contributions  to  the 

j=l 

internal  energy. 


A  frequently  made  approximation  is  (Fumi  and  Tosi,  1962) 


E  =  Eg(V)  +  Ey(T)  (10) 

and 

E  -  Eq(V)  +  Et(T)  (11) 

This  division  of  total  internal  energy  in  a  volume -dependent  term  and  a 
temperature -dependent  term  is  called  Hildebrand's  approximation. 


In  this  report  the  contribution  of  electrons  to  the  internal  energy  will 
be  largely  ignored.  It  should  be  mentioned,  however,  that  for  metals  at  low 
temperatures,  electron  heat  capacities  predominate  over  the  lattice  heat 
capacities  (Gopal,  1966,  p.  62). 

For  a  solid  the  Helmholtz  free  energy  function  (Fumi  and  Tosi,  1962) 

F  =  E  -  TS 


may  be  written 


F  =  EQ  +  ET  -  TS 

and  F^,,  the  thermal  free  energy,  defined  as 

Ft  =  Et  -  TS 


(12) 
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Thus,  the  pressure  at  any  given  temperature  becomes 


or 

P(V ,  T)  =  PQ(V)  +  Pt(V,T)  (13) 


where  P^,  may  be  defined  as  the  thermal  pressure  which  becomes  zero  at 
absolute  zero.  Equation  12  may  be  written  in  terms  of  vibrational  energies, 
1*6«| 


and 


or 


(14) 


where  Py  may  be  defined  as  vibrational  pressure.  Py  does  not  become 
zero  at  absolute  zero.  The  term  -dEg/dV,  which  it  would  seem  reasonable 
to  designate  as  some  sort  of  pressure,  does  not  appear  to  have  a  name  in 
the  literature.  We  will  call  it  the  static  pressure. 


A  few  more  of  the  features  of  the  zero  degree  isotherm  should  be 
pointed  out:  The  term  E^  in  Figure  1-1  represents  the  energy  that  would  be 
required,  at  absolute  zero,  to  separate  the  atoms  of  the  solid  to  infinity. 

It  is  called  cohesive  energy  and  is  approximately  equal  to  the  heat  of  evapor¬ 
ation  at  absolute  zero,  (Zel'dovich  and  Rayzer,  1965,  p.  843).  It  should  be  noted 


« 

\ 

«•% 
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that  the  term  cohesive  energy  also  is  used  in  the  literature  for  E^.  In  this 
report,  however,  cohesive  energy  will  be  reserved  for  E^,. 

c 


The  cohesive  energy  can  be  given  by 

00 


L 


podV=  -Ec 


The  fundamental  equation  of  thermodynamics  which  is  applicable  to 
any  reversible  process,  or  irreversible  process  providing  that  the  states 
considered  are  equilibrium  states,  is 

dE  =  TdS  -  Pd  V  (15) 

At  zero  degrees  absolute  T,  dS  and  TdS  are  equal  to  zero,  so  that 


dE  =  -PdV 


(16) 


and 


(dP/dV).  =  (SP/dV), 


(17) 


Hence,  the  zero  degree  isotherm  is  also  an  adiabat  for  which  dS  =  0. 


The  compressibility  at  zero  pressure  and  ordinary  temperatures  is 
approximately  equal  to  the  compressibility  at  zero  pressure  and  absolute 
zero  (Zel'dovich  and  Rayzer,  1965,  p.845),  i.e., 


Thus,  the  speed  of  sound,  considered  as  a  purely  thermodynamic  quantity, 


is  determined  by  the  slope  of  the  PV  zero  degree  isotherm. 
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The  purpose  of  the  preceding  discussion  is  to  bring  together  in  one 
place  the  various  definitions  and  formulations  of  thermodynamic  terms 
pertaining  to  solids  found  scattered  in  the  literature  as  an  aid  for  quick 
reference. 
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Section  2 

THERMODYNAMICS  OF  EQUATION  OF  STATE 

A  thermodynamic  state  is  specified  in  terms  of  the  independent  vari¬ 
ables  of  the  system.  The  number  of  independent  variables  necessary  to  specify 
the  state,  however ,  must  be  decided  by  the  experimenter  (Reiss,  p.  5,  1965). 

For  example,  usually  only  two  state  variables  are  necessary  to  specify  the 
state  of  a  solid,  but  if  the  surface  area  becomes  appreciable,  as  in  a  powder, 
another  state  variable  is  needed. 

A  system  is  completely  defined  thermodynamically  if  the  entropy  is 
known  as  a  function  of  the  extensive  parameters.  For  a  closed  system  in  which 
the  amount  of  single -phase  matter  remains  constant  and  which  can  be  specified 
as  a  function  of  two  extensive  parameters,  the  entropy  can  be  written  as 

S  =  S(E,  V)  (1) 

Equation  1  may  be  written  in  the  alternate  form, 

E  =  E(S,  V)  (2) 

The  differential  form  of  Equation  (2)  determines  changes  of  state: 

dE  =  (dE/aS)ydS  +  (0E/dV)gdV  (3) 

The  partial- derivatives  in  the  above  equation  have  the  following  meanings 
(dE/dS)y  =  T  and  (aE/aV)g  ■  -P 
The  above  partials  were  derived  from 

dE  =  TdS-PdV  (4) 
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which  is  called  the  fundamental  equal  ion.  Equivalent  expressions  of  Equation 
(4)  are: 


dH  =  VdP  +  TdS 

(5) 

dF  =  -PdV  -  SdT 

(6) 

dA  =  VdP  -  SdT 

(7) 

Because  Equations  { 1 )  and  (2)  contain  all  the  information  necessary  to 
completely  define  a  system  the rmodynamirally,  they  are  called  fundamental 
relationships  or  complete  equations  of  stale.  Other  complete  equations  of 
state,  corresponding  to  Equations  (5)  through  (7)  are 

H  =  H(S,P) 

F  =  F(V,  T) 

A  -  A(P,  T) 

As  is  evident  in  the  above  discussion,  knowledge  or  specification  of 
two  independent  variables,  S  and  V,  or  E  and  V,  depending  on  whether 
Equations  (l)  or  (2)  are  used,  completely  defines  the  system,  providing, 
of  course,  that  E  or  S,  the  dependent  variables,  are  known  as  a  function  of 
E  and  V  or  S  and  V.  The  quantities  S  and  E,  however,  are  not  measurable 
experimentally.  Variables  that  can  be  measured  experimentally  are  P,  V, 
and  T.  Because  P  and  T  are  functions  of  S  and  V,  specification  of  any  two 
of  the  variables  P,  V,  or  T  also  completely  determines  the  thermodynamic 

i 

state  of  the  system. 

Experimentally,  the  functional  relationships  that  are  usually  determined 

are: 

P  =  P(V,  T) 

and 

E  =  E(V,T) 

Neither  one  of  the  above  two  equations,  however,  supply  all  the  informa¬ 
tion  necessary  to  completely  define  the  system  thermodynamically.  To  do  so, 
both  equations  must  be  used.  In  other  words,  although  specification  of  V  and  T 
completely  determine  the  system,  more  than  one  equation  is  necessary  to  obtain 


2-2 


LMSC/HREC  A784868 


tell  the  thermodynamic  information  about  it. 

Because  thermal,  P=P(V,T),  and  caloric,  E=E(V,T),  equations  of 
state  do  not  supply  by  themselves  all  the  necessary  information,  they  are 
called  incomplete. 


The  preceeding  discussion  of  equations  of  states  is  based  on  similar 
discussions  in  Anderson,  Doran,  and  Fabrenbruch,  1965,  pp.  3-5,  and  Callen, 
I960,  p.  25  and  pp.  33-34. 

The  following  discussion  concerns  itself  with  the  derivation  of  some 
incomplete  equations  of  state. 


An  equation  of  state  may  be  derived  from  the  free  energy  because 

P(V ,  T)  =  -(dF/dV) 


(8) 


Using  the  vibrational  form  of  free  energy,  Equation  (8)  becomes 


dE  (V)  /d  F, 
-P  «  — ~ —  + 


dV 


(9) 


Mr'  “ 


Performing  the  following  manipulations  on  the  term 

easily  shown  that  Equation  (9)  reduces  to  an  equation  similar  to  the 
Mie-Gruneisen  equation  of  state. 


can  be 


The  manipulations 


are: 
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The  thermal  formulation  of  free  energy  lends  to 


P  + 


dE0(V) 

dV 


T 


or 


P  -P0(V) 


(10) 


(H) 


The  terms  7^  and  are  generalized  Griineisen  parameters  equrl  to 
(Fumi  and  Tosi,  1962) 
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Let  us  now  see  what  sort  of  an  equation  results  from  manipulation  of 
the  following  equatidn: 

(ff)  =  -P(S.V)  (14) 

The  manipulations  are: 
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So  Equation  (15)  becomes  (Fowles,  1964,  p.  51 ) 


P  -  P 


y3  /d  e  \  /8Et\ 

>(v)  =  v"  Tywjv  ' 


(16) 


Other  equivalent  ways  of  writing  Equation  (16)  are : 

Y, 


where 


and 


where 


and 


p  '  po  =  T  [E  Eo]  +  x(v- T) 


X(V,0) 


,  /9X\  /E-Eo\/8y3\ 

’  \0T/V  ■  y  v  )\ST/^ 


P  =  -y  E  +  R(V ,  T) 


R(V,  0)  =  - 


dE0  r3(v,o) 
_  _ 


/8R\  E  (d  y3\ 

\Ff/v  -  ■  v 


(17) 


(18) 


Equations  (16),  (17),  and  (18)  are  general  Mie -Gruneisen  equations  of 
state.  In  these  general  Mie-Gruneisen  equations  Y^  is  a  function  of  two 
state  variables,  usually  given  as  V  and  T. 
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It  can  also  be  expressed  in  the  following  ways:  (Anderson,  Doran,  and 
Fabrenbruch,  1965,  p.  6) 


y3(x,  Y) 


(9P/9T)v  V(8S/8V)t  ya 

V  (8E/8T)V  =  T(8S/8T)V  =  £~C^, 


The  terms  X  and  Y  in  the  above  expression  stand  for  the  variables  Y ^  is  a 
function  of.  For  example, 


y3(E,  v) 


To  sum  up  the  discussion  concerning  the  derivation  of  incomplete 
equations  of  state,  Equations  (10)  and  (11)  which  involve  y^  and  Y^  are 
equations  of  isotherms  while  Equations  (16),  (17)  and  (18)  which  involve  Y ^ 
are  equations  of  isentropes. 
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Section  3 


A.  THE  MIE -GRUNEISEN  EQUATION  OF  STATE  AND  THE  GRUNEISEN 
CONSTANT 


The  assumption  that  Y  is  a  function  of  volume  only  restricts  the 
fundamental  dependencies  of  T,  S,  and  E,  and  makes  possible  the  derivation 
of  a  useful  equation  of  state.  The  following  discussion  should  make  obvious 
these  points. 


If  Xj  is  a  function  of  volume  only, 


then  integration  of  Xj  = 


yields 


f  dT 

J  T 

V3  ^  +  in  0(S) 

In  T 

X3  ^  +  in  *(S) 

fy  £¥ 

T 

m 

-%e" 

'  3  v  =  «V> 

T(S,  V) 

=  0(S)  0(V) 

s  ■  4m) 


The  form  of  E  is  obtained  by  integration  of 


or 


T  *(H)V  * 

E  -  E0<V)  =  0(V)  fm  dS 


(1) 
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Now 


=  P 


so  differentiation  of  Equation  (1)  leads  to 


or 


P  + 


P  -  PQ(V) 


dfl(V) 

dV 


J  <t>(S)  dS 


0 


yG(V) 

V 


s 

4>(S)  dS 


0 


where  7L 
u 


is  Griineisen's  constant  aid  is  defined 


by 


yGr/) 


d  In  ft(V) 
din  V 


(2) 


(3) 


Equation  (2)  finally  can  be  reduced  to 


yG(v) 


P  -  P0(V)  *  -if—  (E  -  E0(V)) 


or 


yr(v) 

p  -  p0(v)  =  et 


Utilising  the  vibrational  formulation  of  energy, 


E  =  ES(V)  +  EV(V,T) 


leads  to 


dE, 


P  + 


dV 
d  E, 


P  + 


yr(v) 

-V<E-v 

_  yG<V>  - 


dV 


(4) 

(5) 

(6) 

(7) 


Equations  (4)  and  (5)  or  (6)and(7)  are  the  Mie-Griineisen  equations 
of  state. 


The  preceding  discussion  was  taken  mainly  from  Anderson,  Doran,  and 
Fabrenbruch,  1965,  pp.  5-6,  and  Fowles,  1964,  pp.  49-53. 
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The  preceding  derivation  of  the  Mie -Gruneisen  equation  is  not  the 
original  derivation  of  Gruneisen.  Originally,  Gruneisen  derived  Equation  (7) 
from  force-field  theory  (Partington,  1957,  p.  346).  Various  other  derivations 
from  Debye's  theory  and  from  entropy  have  also  been  presented.  It  would  not 
be  worthwhile,  for  the  purposes  of  this  report,  to  go  into  the  details  of  these 
derivations.  However,  a  point  of  definition  of  y ^  which  appears  in  these 
derivations  is  important. 

In  statistical  mechanics  the  isothermal  equation  of  state  may  he  repre¬ 
sented  as  (Slater,  1939,  p.  21 7) 


P  -  P 


hv. 

_ 1 

h  v. 


I 


where  v.  represents  the  vibrational  frequency  of  the  j  normal  mode  and 

dtnv  . 


/ornv.i 

-  ■  (aKv). 


Gruneisen  assumed  that  all  of  the  were  equal  to  each  other  and  to 
a  constant  y  (Slater,  1939,  p.  219).  Under  this  assumption 


S 

A. 


i 

Am 


(d in  v ,  v 
dtoiV  ) 


d  fn  v . 

_ 1  = 

d  fn  V 


y3  =  *G 


According  to  the  Debye  theory,  the  limiting  frequency  vmax  defines  the 
frequency  spectrum.  The  frequencies  of  all  the  oscillations  change  in  the 
same  ratio  as  the  change  in  vmax*  (Slater,  1939,  p.  238) 


Thus 


dfnv 


r  = 


max 


d  fn  V 
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In  practical  applications  of  the  Mie-Gruneisen  equation,  the  Gruneisen 
parameter  y^  has  been  used  rather  indiscriminately.  Fumi  and  Tosi,  1962, 
have  noted  the  conditions  under  which  its  use  is  valid.  The  following  is  a 
discussion,  taken  from  their  paper,  of  the  use  of  y^  in  equations  of  state. 


Because  Y^  is  usually  the  easiest  to  determine  experimentally,  it  has 
often  been  substituted  for  Y j  and  Y^.  Such  a  substitution,  however,  is  valid 
only  if  or  Y ^  do  not  change  with  temperature,  as  seen  from  the  following 
equations  which  relate  Xj  with  ^  and  Y^: 


It  does  not  follow  that  because  Y ^  does  not  change  with  temperature 
that  and  Y ^  also  do  not.  It  can  be  shown,  however,  that  if,  in  a  tempera¬ 
ture  range  from  0  K  upwards, 


as  would  be  the  case  if  Xj  were  a  function  of  volume  only,  then  y^  =  Y ^  in  the 
same  temperature  range.  This  equality  also  applies  if,  from  0°K  upwards, 

CV(V.T)  =  Cy^) 
and  the  entropy  at  0°K  is  a  constant. 


In  other  words,  to  summarize  the  preceding  discussion,  if  Y ^  or  Y ^  are 
functions  of  volume  only,  then 


or 


3-4 


LMSC/HREC  A784868 


Because  Y^  is  a  function  of  volume  only  in  a  certain  range  of  temperature 
does  not  mean,  however,  that  Equations  (4)  and  (5)  or  (6)  and  (7)  are  necessarily 
valid.  Equations  (4)  and  (5)  are  valid  if  Y 3  is  a  function  of  volume  only  in  the 
temperature  range  from  zero  degrees  absolute  upwards,  or  if 

Cy(V,  T)  =  Cy^TV))  (0°K  upwards) 

in  the  same  temperature  range  and  the  entropy  at  zero  degrees  absolute 
is  constant. 


If  the  energy  can  be  expressed  according  to  Hildebrand's  approximation: 


E  =  Eq(V)  +  Et(T) 
or 

E  =  Eg(V)  +  Ey(T) 


then  the  relationship  between  the  various  Ys  are 


and 


^(V.T)  = 

TCy  y3(V,T) 

(10) 

r2(v,T)  = 

TCy  y3(V,T) 

(11) 
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To  convert  Bridgman's  data  to  Staler's  notation  (Overton,  1962,  p.  118): 


V0T  =  V00(I+5V 


a  =  a 


l/‘+a0 


b  =  V  ta0 


Equation  (1)  can  be  rewritten  in  the  following  form 


P  =  P0(T)  +  P,(T)  [jv0()  -  V)/V00[  +  P2(T)  |<V00  -  V). 


»]  +  P2<T>  [<V00  -  V»Oo]' 


(6) 


where  Pq(T)  becomes  zero  at  absolute  zero.  In  terms  of  Equation  (6)  the 
heat  capacity  at  constant  volume  is  (see  Slater,  1939,  p.  204) 


cv  -  cv  V00T 


[d2Px/dT2  [<V00  -  V)/V00j  +  1/2  d2p/dT2  [(V0#-V)/V00]2 


+  l/3d2P2/dT2[(V00.V)/v00j3... 


(7) 


O  . 


where  Cy  is  the  heat  capacity  at  volume 


00- 


b.  Gruneisen  Parameter  y^ 

As  mentioned  previously,  the  Gruneisen  parameter, y^,  is  usually 
determined  experimentally  from  the  relationship 


y3(v.T)  =  av//9Tcv 


(8) 


4> 

This  Pq(T)  is  not  the  same  as  Pq(V)  used  previously.  The  latter  represents 

the  elastic  pressure  on  the  eero-degree  isotherm  while  the  former  is  the 
pressure  required  to  reduce  the  volume  to  VQ0.  To  avoid  confusion  PQ(T) 
will  be  called  P^(T)  in  this  paper. 
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or  alternatively 

y3(V,T)  =  CtV/0 sCp  (9) 

where  0C  is  the  adiabatic  compressibility 

d 


{  =  -(l/V)(dV/dP)s,  fisCp  =  flfCv) 


Like  any  other  state  variable,  the  Gruneisen  parameter  is  a  function 
of  two  independent  variables.  If  one  independent  variable  is  set,  then  the 
parameter  becomes  a  function  of  the  other.  Because  the  independent  variables 
most  convenient  to  work  with  arc  pressure  and  temperature,  Equation  (8)  could 
just  as  well  be  written  in  terms  of  P  and  T.  Such  an  equation  can  be  written 


y3(P.T) 


q(P,T)  V(P,  T) 
fiiP.T)  Cy(P,  T) 


(10) 


If  P  is  set  at  aero,  Equation  (10)  becomes  &  function  of  temperature  alone. 
Equation  (10)  thus  gives  us  the  temperature  dependency  of  y3  along  an  isobar. 
Because  the  volume  expansion  of  a  solid  at  zero  pressure  is  relatively  slight  - 
aluminum  increases  its  volume  by  only  some  5%  from  20°C  to  600°C  -  Equation 
(10)  is  only  slightly  volume -dependent. 


In  principle,  the  volume  dependency  of  y^  could  be  checked  by  holding 
temperature  constant  and  letting  volume  or  pressure  vary.  This  procedure 
would  require  thermal  expansion,  heat  capacity,  and  compressibility  data  at 
various  pressures.  Unfortunately,  such  data  is  not  generally  available. 


Barron,  Leadbetter,  and  Morrison,  1964,  p.  75,  have  made  a  rough  esti¬ 
mate  of  the  volume  dependency  of  y^  from  Equation  (8)  for  NaCl  and  KC1. 
Their  thermodynamic  equation,  derived  from  Equation  (8),  for  this  depend¬ 
ency  is 


(dtn  y3 

)  _L  U 

1  p<*rv>l  T 

d(Vcr) 

J_| 

[•*)  j 

(  InV 

>T  *r  r 

y  L  Jp  cp 

l9T  Jp 
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At  T  300°K  they  found  (pressure  presumably  zero) 


( 


a/ny3^ 
dlnV  )T 


0.3  (NaCl) 
1.9  (KC1) 


It  is  pointed  out,  however,  that  the  accuracy  of  the  above  results  may 
be  in  error  by  +  100%,  or  more,  for  NaCl  and  +50%  for  KC1. 

Making  various  approximations  and  interpolations  of  experimental  data, 
Ku  (1967)  used  Equation  (10)  to  calculate  Y ^  values  at  different  temperatures 
and  pressures  for  a  number  of  different  plastics.  His  results  are  given  in 
Tables  4-1  through  4-4. 


c.  Determination  of  Gruneisen  Parameter  from  Thermal  Expansion  and 
Heat  Capacity  Data 


Although  heat  capacity  and  thermal  expansion  data  at  zero  pressure  are 
readily  available  for  most  elements  and  compounds,  compressibility  data 
over  extensive  temperature  ranges  is  very  limited.  Thus,  Equation  (10) 
cannot  be  used  in  most  cases  to  test  the  temperature  dependence  of  y y 
Another  equation,  however,  was  derived  by  Gruneisen  which  does  not  require 
compressibility  data  over  an  extended  temperature  range.  The  derivation  of 
this  equation  is  based  on  the  expansion  of  the  term  V  dEc/dV  in  the  Mie- 

U 

Gruneisen  equation.  At  zero  pressure  (Partington,  1962,  p.  351) 

dES 

-V-d VT-  ■  VCEV 


and  after  expansion  of  V  dEg/dV  in  powers  of  (V  -  Vqq) 


1  +  (n  +  m)/3 


*GEV 


(12) 
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GRUNEISEN  PARAMETER 
Plastics 

Polyethylene  -  low  density 

Polyethylene  -  high  density 

Teflon 

Nylon 

Kel-F 

PMMA 

Polyester  (Sectron  5003)° 
Polystyrene 

Polypropylene 

Cellulose  acetate  Butyrate 

Ethyl  Cellulose 

Polyvinyl 
Chloride  -  Pure 

-  10%  Plasticizer 

-  20%  Plasticizer 

-  30%  Plasticizer 
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Table  4-1 

OF  SOME  COMMON  PLASTICS  (Ku,  1967) 


Gruneisen 

Parameter 

Pressure 

atm 

Temp 

°c 

1.213 

10,000 

20 

0.812 

2,000 

20 

0.941 

10,000 

20 

0.796 

1 

20 

0.305 

10,000 

20 

0.285 

1,000 

21.1 

1.04 

10,000 

20 

1.44 

1 

— 

0.995 

1,910 

29 

1.80 

618 

23 

0.595 

1,910 

181 

0.480 

1,910 

196 

0.646 

1,000 

0 

0.865 

1,000 

20 

0.819 

1,000 

■i : 

0.866 

1,000 

HI 

1.32 

1,000 

0.910 

1,000 

1.16 

1,000 

0.913 

1,000 

H 
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Table  4-2 

GRUNEISEN  PARAMETER  y,  OF  TEFLON  (Ku,  1967) 


Pressure,  atm 


Temperature  °C 

20 

60 

80 

.935 

.768 

.374 

.290 

.651 

.698 

.411 

.355 

.591 

.611 

.469 

.417 

.547 

.577 

.505 

.459 

.568 

.520 

.538 

.505 

.598 

.557 

.561 

.533 

.615 

.576 

.577 

.556 

.656 

.627 

.611 

.581 

.724 

.680 

.645 

.611 

.819 

.774 

.720 

.657 

.941 

.866 

.808 

.741 

4,000 

5,000 

6,000 

7,000 

8,000 


10,000 


Table  4-3 

GRUNEISEN  PARAMETER  X,  OF  POLYETHYLENE  (Ku,  1967) 


m _ _ _ a. _ 


Pressure,  atm 


20 


.177 
.990 
.664 
.657 
1 .654 
1.559 
1.597 
1.398 
1.332 
1.216 
1.213 


Temperature  C 


40  60 


.815  .539 

.431  1.041 

.164  .885 

1.172  .927 

1.226  .977 

1.198  .959 

1.254  1.015 

1.121  .920 

1.088  .896 

1.019  .840 

1.026  .846 


1.047 

1.172 

1.074 

1.271 

1.330 

1.269 

1.393 

.776 
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where  n  and  m  are  constants  associated  with  the  repulsive  and  attractive 
forces.  Equation  (12),  after  some  simplifying  assumptions,  can  be  written 
(Gruneisen,  1926,  p.  46) 

V. 


OT 


00 


00 


<VG+f)  Ev 


(13) 


where  Q  = 
o 


00 


y  c.  fit 


Cp00 


Differentiation  of  Equation  (13)  with  respect  to  temperature,  assuming 


As  pointed  out  by  Hume-Rothery,  1945,  p.  210,  the  term  (y^. 
increasingly  important  at  high  temperatures. 


2  EV 

+  -j)  q —  becomes 
Uo 


To  a  first  approximation  Equation  (14)  can  be  written 

vooa 


*00  cv 


(15) 


The  assumed  temperature  independence  of  y^-.  can  thus  be  tested  by 
using  O  and  Cy  values  at  various  temperatures  and  rero  pressure. 


d.  Slater  and  Dugdale -MacDonald  Relations 


Because  the  volume  dependencies  of  y^.  and  y^  are  of  great  interest, 
and  the  data  for  their  calculation  from  thermodynamic  equations  are  in  short 
supply,  great  use  has  been  made  of  the  Slater  relationship  (Slater,  1939,  p.  239): 

ys(V)  =  -  2/3  -  1/2  V  (d2P/aV2)r/tdP/aV)T  (16) 
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This  equation  was  derived  hy  assuming  that 

VS<V>  =  -  d  in  vmax/d  In  V 

where  v"rnax  *•  Debye's  limiting  frequency.  It  is  also  assumed  that  the 
Poisson  ratio  is  independent  of  volume. 

The  Gruneisen  parameter  determined  by  Equation  (16)  is  labeled  y _ 

o 

because,  as  Collins,  1962,  p.  323,  points  out,  y g  can  be  identified  with  y^. 
obtained  from  thermal  expansion  only  if  Poisson's  radio  is  independent  of 
volume.  Usually  this  is  not  the  case. 


At  absolute  zero  and  zero  pressure,  the  volume  is  Vqq  and  Equation 
(16)  takes  the  form 


V 


s 


(17) 


where  P°  and  P°  represents  the  temperature  dependent  parameters  at  T  =  0°K. 
Because  to  a  good  approximation  P°  and  P^  can  be  taken  as  equal  to  Pj  and 
P^  at  room  temperature  (Kachhava  and  Saxena,  1965,  p.  986)  and  thermal 
volume  expansion  is  assumed  negligible  (Slater,  1940,  p.  745),  Slater's  relation 
can  be  written 


y 


s 


2 

3 


(18) 


or,  in  Bridgman's  notation 


b  2 

yS  "  2*3 

a 

for  room  temperature  and  zero  pressure. 


(19) 


Gillvarry,  1957,  p.  1258,  using  a  generalized  isothermal  equation  of 

state,  showed  that  even  though  yg  is  temperature  dependent,  this  dependency 

is  not  manifest  at  zero  pressure.  At  zero  pressure  yc  is  strictly  constant 

o 

at  any  temperature.  Thus,  Slater's  relation  apparently  cannot  be  used  to  test 
the  temperature  dependency  of  yg  at  zero  pressure. 
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If  the  preceding  assumptions  are  granted,  Equation  (16)  can  be  used  to 
determine  the  volume  dependence  of  y^.  In  Bridgman's  notation,  Equation 
(16)  becomes 


ys(V) 


v/voo 

[  ‘  +  2  -r  <l  -  v/voo>| 

1  a  ' 


(20) 


To  the  best  of  our  knowledge  the  above  equation  has  not  appeared  as 
such  in  the  literature  nor  has  a  critical  evaluation  of  its  validity. 


There  appears  in  the  literature  another  equation  for  Slater's  relation. 
This  is 


(21) 


Although  no  discussion  has  been  found  in  the  literature  about  the  condi¬ 
tions  under  which  Equations  (16)  nr  (21)  should  be  used,  some  pertinent  points 
have  been  found: 


Slater's  derivation  of  Equation  (16)  is  based  on  the  relationship  (Slater, 
1939,  p.238) 


v 

max 


(22) 


where  v^  and  v^  are  the  velocities  of  the  longitudinal  and  transverse  waves. 
The  wave  velocities  are  related  to  compressibility  by  the  following  equations. 

v  =  J  3<i  and  v  =  JjQgf 

v|  \  1  -<r  )  and  Vt  y  2fip(  1  +cr) 

where  a  is  Poisson's  ratio. 
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Dugdale  and  MacDonald,  1953,  p.  834,  state  that  the  propagation  of  a  sound 
wave  (presumably  under  zero  applied  pressure)  at  low  temperatures  will  be 
isothermal,  while  at  higher  temperatures,  including  room  temperature,  the 
process  will  be  practically  adiabatic.  Gilvarry,  1956,  p.  333,  on  the  other 
hand,  points  out  that  at  low  pressure  the  adiabatic  or  isothermal  compres¬ 
sibilities  are  for  practical  purposes,  identical. 


It  is  apparent  that  more  study  of  this  matter  is  required  to  resolve 
the  confusion.  Of  course,  at  zero  degrees  absolute,  no  ambiguity  exists 
because,  at  T  =  0°K, 


Dugdale  and  MacDonald  (1953)  objected  to  Slater'B  relation  because 
they  claimed  that  the  effect  of  finite  strain  was  neglected.  They  proposed 
the  following  relationship  instead 


v  az(P  v2/3)/av2  .  l 
2  d(P  v  2/3)/dV  3 


(23) 


Their  reasoning  was  shown  to  be  in  error  by  Gilvarry  (1956)  and  Barron  (1957). 
However,  Equation  (23)  has  shown  good  agreement  for  zero  pressure  tests 
(Rice,  McQueen,  and  Walsh,  1958,  p.  44).  Generally,  there  is  no  firm  basis 
for  preferring  Equation  (23)  over  Equation  (16)  (Doran,  I960,  p.  29). 


e.  Gruneisen  Parameter  y ^ 

If  a  means  is  available  for  obtaining  P^,(V,  T),  then  y £(V ,  T)  can  be 
determined  from 

y,(v,  T) 

P  =  — -  E 

T  V  T 


The  term  P^,  can  be  easily  obtained  if  the  P-V  zero  degree  isotherm 
and  P-V  isotherms  are  available.  From  such  data,  P^  is  obtained  by 
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P^,  =  P  -  Pq.  Unfortunately,  compressibility  data  is  in  short  supply,  so 
some  other  expedient  for  obtaining  P^,  is  desirable. 

Bee  croft  and  Swenson  (1961)  found  that  if  the  isothermal  compressibi¬ 
lity  of  a  substance  is  linear  with  volume  and  independent  of  temperature,  the 
determination  of  P^,  is  greatly  facilitated.  The  following  will  briefly  review 
the  approach  of  Beecroft  and  Swenson. 

If  the  compressibility  is  a  function  of  volume  alone  and  independent  of 
temperature,  then  P^.  becomes  a  function  of  temperature  alone.  This  becomes 
obvious  if  one  integrates  the  expression  for  compressibility,  i.e., 


(24) 


Thus,  if  the  compressibility  is  a  function  of  volume  alone,  Po(V)  and  P^.(T) 
can  be  evaluated. 


Furthermore,  if  P^,  is  a  function  of  temperature  alone  then 

Id P\  a  dPT 
Wly  "  P  “  dT 

so  that  a/p  is  a  function  of  temperature  alone.  Thus  a/P  found  at  P  =  0  can 
be  considered  as  a  constant  for  a  given  temperature. 

The  expression  for  E^,  is  found  by  integrating 
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Thus, 


(25) 


and 


Subtracting  Equations  (25)  and  (26)  gives 


(26) 


Et  =  E(V,  T)  -  Eq(V)  =  E(Vqt,T) 


eoo  +  t|(vpt-vot> 


(27) 


The  terms 

e(vot’  t)  ■  Eoo 

can  be  replaced  by  H  at  P  =  0  and  T.  This  follows  from 


at  P  =  0 


and 
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The  term  Eqq  is  arbitrarily  set  equal  to  zero.  The  last  two  integrals 
of  Equation  (27)  can  be  evaluated  from  experimental,  compressibility  data. 

Beecroft  and  Swenson  found  for  sodium,  that  the  compressibility  could 
be  assumed,  to  a  first  approximation,  as  independent  of  temperature  and  a 
linear  function  of  volume.  Their  experimental  data  is  shown  in  Figure  4-1. 


Figure  4-1  -  The  Isothermal  Compressibility  of  Sodium  as  a  Function  of  Molar 
Volume 


From  this  data  and  heat  capacity  data  at  zero  pressure  the  temperature 
and  volume  dependence  of  was  calculated.  The  results  are  shown  in 
Figure  4-2. 
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Variation  of  the  Gruneisen  Constant 


Figure  4-2  -  The  Volume  Variation  of  the  Gruneisen  Constant  as  Given 
by  Various  Definitions 

B.  UNDER  DYNAMIC  CONDITIONS 

a.  Method  of  Rice,  Walsh,  and  McQueen 

In  the  shock  wave  determination  of  the  Gruneisen  constant,  the  form 
of  the  Mie-Gruneisen  equation  of  state  usually  assumed  is  the  thermal 
formulation 

P  -  PQ  =  <yG/V)(E  -  Eq)  (28) 

Equation  (28)  combined  with  the  shock  Hugoniot, 

eh  -  eoh  *  <V2«voh  -  v>  <2’> 

where  E^  and  are  conditions  ahead  of  the  advancing  shock  and  P^  the 
pressure  behind  the  shock,  gives 
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*0  ■  V<PH 


-  P 


o)/[<ph/2)(voh 


-  V)  +  E 


OH 


(30) 


where  Eq  = 


dV  (E 


00 


the  energy  at  zero  degrees  absolute  and  zero 


pressure  is  arbitrarily  assumed  zero).  The  term  E^  is  found  by  integrating 
the  Debye  specific  heat  curve  from  zero  degrees  absolute  to  the  temperature 
at  the  initial  shock  conditions. 


It  is  further  assumed  that 

yG  =  -  1/2  [(d(P0V2/3)/dV2)(d(P0V2/3)/dV)]  -  1/3  (31) 

This  latter  equation  is  the  Dugdale -MacDonald  modification  of  Slater's 
formula  assumed  to  be  valid  along  the  zero  degree  isotherm.  Thus, 

-  V/2  [(d(P0V2/3)/dV2)/(dP0V2/3/<iV)]  -1/3 

*  V<PH  -  P0>/  [(PH/2)(V0H  -  V>  +  E0H  '  Efl]  <32> 

Knowing  the  initial  conditions  of  specific  volume,  compressibility  at 
zero  degrees  absolute  and  zero  pressure,  and  specific  internal  energy  E^^ 
(relative  to  Eqq),  Equation  (32)  can  be  integrated.  The  result  is  the  zero 
degree  isotherm.  From  this  zero  degree  isotherm,  y^  is  evaluated  using 
Equation  (31 ). 

The  preceding  discussion  was  taken  mainly  from  Rice,  McQueen,  and 
Walsh,  1958,  and  Deal,  1962. 

Rice,  et  al,  1958,  p.  57,  present  their  calculations  of  y^  in  the  form 

*G(V)  =  V0H  +a^+b^2  +  c^3  (33) 

where  H  =  P/°qh  "  A,  B,  and  C  are  constants,  and  is  the  value  of  y^at  zero 


4-16 


LMSC/HREC  A784868 


pressure.  The  values  of  Vq^  were  choosen  to  correspond  to  those  obtained, 
at  aero  pressure,  from  the  equation 


y 


VO 

/?cv 


(34) 


To  check  the  consistency  of  the  preceding  method,  the  Dugdale- 
MacDonald  y  is  found  at  sero  pressure  from  the  relationship 


y  =  2S  -  1  (35) 

where  S  is  the  constant  in  the  linear  relationship  between  shock  velocity  and 
particle  velocity,  i.e., 


=  CA  +  SV, 


Equation  (35)  represents  the  Dugdale -MacDonald  parameter  at  aero  pressure 
as  derived  entirely  from  shock  data.  The  derivation  of  Equation  (35),  which 
is  not  discussed  here,  is  straightforward  (see  Deal,  1962,  p.  206).  The  param¬ 
eter  determined  from  Equation  (35)  may  now  be  compared  to  that  determined 
by  Equation  (34).  Such  a  comparison  shr  ws  that  the  two  y's  agree  fairly  well 
for  many  elements  (McQueen,  1964,  p.  68). 


b.  (dE/8P)y  As  sumed  Constant 


If  the  assumption  is  made  that  ( 6E/dP)y  is  independent  of  pressure  and 
volume  in  the  vicinity  of  the  Hugoniot,  and  y  ^  is  a  function  of  volume  only, 
then  y^.  can  be  determined  from  the  equation  (Deal,  1962,  p.  206,  and  McQueen 
and  Marsh,  I960,  p.  1267): 


rG(v) 


(36) 
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The  term 
data  by 


is  determined  from  aero  ressure, 


the  r  modynamic 


(dE\  _  (aT)v  _  CV 

Mv=(|E)v=  « 

Deal,  1962,  p.  206,  states  that  the  above  procedure  is  justified  because 
it  reproduces  well  equations  of  state  obtained  by  more  complex  procedures. 
A  comparison  of  the  results  obtained  by  this  method  with  other  methods, 
however,  has  not  been  found. 

C.  y from  Speed  of  Sound  Measurements  behind  Shock  Waves 

An  adiabatic  equation  of  state  may  be  written  in  terms  of  the  Hugoniot, 
just  as  an  isotherm  was  in  Section  4  B.  a.  Consider  Figure  4-3. 


Figure  4-3  -  P  vg.  V  Curve  to  Illustrate  Relationship  Between  Sound  Speed 
Behind  a  Shock  Wave  and  Intersection  Between  Adiabats  and 
Hugoniots 
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At  volume  the  pressure  on  the  Hugoniot  can  be  given  by 


PH  ~  P0  =  vf  <EH  -  E0> 


while  that  on  the  isentropc  by 


Thus, 


Or,  generally 


PA  ■  po  *  T  <ea  -  Eo» 


PA  "  PH  V?  (EA  -  EH> 


PA  -  PH  ‘  V  (EA  *  EH> 


Differentiation  of  the  above  equation  gives  (Deal,  1962,  p.  205) 
dEH 

PA  =  '  TV“  +  (PH  '  PAJ 


m.  v 

rdpH 

dPA' 

LyJ  y 

dV 

dV 

m 

(37) 


At  volume  V ^  where  the  Hugoniot  and  adiabat  intersect,  Equation  (37) 
becomes 


P  =  P 
H  A 


dV  y  dV  dV 


or 

dPH 

dpA] 

v 

_  .  Lav 

dV, 

¥ 

d£ 


P  T  Y  + 


H 


(38) 


H  1  dV 


dE 

To  find  ■  jy  ■  ,  the  equation  of  the  Hugoniot  is  differentiated,  i.e. 
eh  •  E0H  ■  T1  <voh  •  v> 
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Ifa  .  i 

dV  ‘  2 


dV 


So  that  Equation  (38)  becomes,  at  volume  (McQueen,  1964,  p.  70) 


vv) 


dP 


H 


dV 


dP/ 

dV 


2V 


dP 


PH  + 


H 


dV  <V0H-V> 


(39) 


The  speed  of  sound,  it  may  be  recalled,  is  given  by 


Thus,  if  the  speeds  of  sound  along  a  Hugoniot  can  be  measured  and  the 
slopes  of  the  Hugoniot  as  functions  of  volume  known,  y^fV)  can  be  calculated 
with  the  aid  of  Equation  (39). 


The  preceding  discussion,  taken  from  McQueen,  1964,  p.  70,  and  Deal, 
1962,  p.  205  represents  the  general  basis  of  the  method.  The  actual  deter¬ 
mination  of  y  from  sound  speeds  behind  shock  waves,  however,  (Al'tshuler, 
Kormer,  Brazhnik,  Vladimirov,  Speranskaya,  and  Funtikov,  I960;  see  also, 
Anderson,  1965)was  approached  in  a  slightly  different  manner  apparently, 
because  Hugoniot  slopes  obtained  from  dynamic  measurements  are  inaccurate. 
In  the  experimental  paper  (Al'tshuler  et  al,  1960),  it  is  reasoned  that  since 


and 


then 


or 


H 


=  +  P, 


>  0 
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where  is  the  isentropic  sound  speed  along  the  Hugoniot  and  Cq  the  sound 
speed  at  absolute  zero.  The  sound  speed  Cq  may  be  written 

dPn 

(40) 


d  a 


where  PL  is  the  density  of  the  initial  Hugoniot  state  and  a  -  P/Pn-  If  expression(40) 
0  2  2  U 

is  integrated  using  instead  of  Cq,  i.e., 


Po*®' 


■f 

500 


P0  CH  d<T  = 


I«t)  - 


KcToq)  =  I (c,  * 


00' 


where  0qq  =  P/Pqq>  Pqq  being  the  density  at  T  =  0°K  and  P  =  0,  then  the  function 
I(<r,  <7qq)  provides  an  upper  limit  to  the  zero  degree  isotherm. 


The  Gruneisen  parameter  is  then  determined  by  using  a  Mie -Gruneisen 
equation  for  P^  which  includes  electronic  contributions.  ForPQ,  I(a,  <Tqq)  is 
substituted. 


The  Gruneisen  parameters  obtained  in  this  manner  are  20  to  25%  lower 
that  those  obtained  by  the  method  of  Rice,  Walsh  and  McQueen  in  which  the 
electronic  contributions  are  taken  into  account. 

d.  Other  Methods 


Other  methods  of  determining  the  Gruneisen  parameter  which  are 
mentioned  in  the  literature  include  varying  the  initial  state  before 
shocking  by  heating  or  by  the  method  of  reflected  shocks.  However,  very 
little  work  appears  to  have  been  done  in  this  area.  These  methods  have  not 
been  investigated  thoroughly. 
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Section  5 


CALCULATIONS  OF  VARIOUS  GRUNEISEN  PARAMETERS  FOR 

MATERIALS  OF  STUDY 


A.  GRUNEISEN'S  PARAMETER  y ^  FROM  HYDROSTATIC  DATA  AT 
ZERO  PRESSURE 


a.  Calculations 


As  discussed  in  Section  4A.b.  ,  the  relationship 


y3(v,  T) 


gy 

PT  Cy 


a  V 


is  usually  used  to  test  the  constancy,  at  aero  pressure,  of  y^  with  tempera¬ 
ture  (the  slight  volume  dependency  of  y3  because  of  thermal  expansion  is 
usually  ignored). 


To  use  the  above  formula,  thermal  expansion,  compressibility,  and  heat 
capacity  data  at  various  temperatures  are  required.  Although  thermal  expan¬ 
sion  and  heat  capacity  data  at  zero  pressure  are  readily  available,  compres¬ 
sibility  data  is  scarce.  For  the  materials  of  this  study,  compressibility  data 
over  an  extended  temperature  range  was  found  only  for  aluminum.  For  the 
other  materials,  with  the  exception  of  the  plastics,  enough  data  was  found  to 
enable  calculations  of  y^  at  room  temperature. 

Because  the  plastics  of  the  study  are  composite  materials  that  have 
varying  compositions  and  porosities,  they  presented  special  problems  which, 
unfortunately,  have  not  yet  been  resolved. 


this  section  the  references  are  numbered. 
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The  results  of  the  calculations  of  y^  from  hydrostatic  data  at  zero 
pressure  are  presented  in  Table  5-1. 


b.  Discussion  of  Results 


Because  of  the  approximations  made  in  the  calculations  and  uncertainties 
in  some  of  the  experimental  data,  it  is  hard  to  state  with  certainty  the  errors 
associated  with  the  calculated  values.  Gschneider  (Reference  12,  p.  416)  states 
that  most  of  the  Gruneisen  parameter  data  agree  with  each  other  only  to  about 
10%. 


What  would  be  desirable  to  increase  the  precision  of  the  above  y^  calcu¬ 
lations  are  more  precise  coefficients  of  thermal  volume  expansion. 


Some  of  the  other  uncertainties  which  influence  the  accuracy  of  the 
above  calculation  are: 

1.  The  data  on  thermal  expansions  and  heat  capacities  of  metals 
have  been  evaluated  for  accuracy.  In  Reference  8  the  data  of 
various  investigators  is  reduced  to  curves  through  "most 
probable  values."  No  such  evaluation  exists,  as  far  as  we 
know,  for  compressibility  data.  The  data  available  on  com¬ 
pressibility  is  limited  and  as  yet  has  not  been  critically 
evaluated  in  a  comprehensive  way.  A  brief  discussion  of  the 
errors  involved  in  compressibility  data  is  given  by  Gschneider 
(Reference  12,  pp.  298-305). 

2.  Although  accurate  data  exist  on  heat  capacities  at  constant 
pressure  of  metals  over  a  wind  temperature  range,  uncer¬ 
tainties  creep  in  when  this  data  is  converted  to  heat  capacity 
at  constant  volume. 

The  Nernst-Lindemann  relationship  below  is  used  for  the 
conversion 


Cp  -  cv  *  A  T  CP 

For  copper,  A  varies  only  by  about  8%  over  a  1200°K  tempera¬ 
ture  range  (Reference  25,  p.  264);  its  variation  for  other  materials 
would  need  to  be  looked  into  in  further  detail. 
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3.  Possible  errors  arising  from  crystal  anisotropy  were  not 
investigated.  It  is  assumed  that  the  data  used  to  calculate 
y j  are  representative  of  isotropic  polycrystalline  material. 


From  Table  5-1,  it  can  be  seen  that  for  aluminum  shows  no  consistent 
variation  with  temperature  from  room  temperature  to  77  3°K. 


Although  the  present  aluminum  values  of  y  are  higher  than  those  given 
by  Moelwyn -Hughes  (Reference  17),  very  little  variation  of  7'  with  tempera¬ 
ture  is  indicated  in  both  sets  of  values.  The  reason  for  our  y  values  being 
higher  than  those  of  Moelwyn -Hughes  lies  in  different  values  of  the  thermo¬ 
dynamic  functions.  However,  our  value  at  293°K,  2.25,  compares  well  with 
the  value,  2.  17,  given  by  Griineisen  in  his  paper  "The  State  of  a  Solid  Body" 
(Reference  11). 

Discussion  of  the  temperature  variation  of  7'  found  in  the  literature 
indicates  that  the  variation  is  greatest  at  temperatures  below  the  Debye 
temperature.  Bijl  and  Pullan  present  the  figure  (Figure  5-1)  on  the  follow¬ 
ing  page  for  the  y  variation  at  low  temperatures  (Reference  3). 

In  contrast  to  the  results  of  Bijl  and  Pullan,  the  experimental  results 
of  Figins,  James,  and  Riley  (Reference  6)  indicate  that  y ^  increases  at  very 
low  temperatures.  Apparently,  the  variation  of  7  ^  with  temperature  is  still 
an  open  question. 


Alder  (Reference  1)  makes  the  statement  that  7^  might  be  expected 
to  decrease  from  about  two  to  its  classical  perfect-gas  value  of  2/3, 
if  the  temperature  is  raised  holding  the  density  constant.  On  the  other 
hand,  Cy  and,  hence,  7'  might  be  nearly  independent  of  temperature. 
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I. 

1. 


■>— —  calculated  with  c<m8*’..t  volume  and  compressibility 
-  — corrected  for  variation  of  volume 
*  experimental  value  given  by  Gruneisen 

a  389°K  for  aluminum 


Figure  5-1  -  Variation  of  Gruneisen  Parameter  at  Low  Temperatures 
for  Aluminum  (Reference  3) 
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The  calculated  y^  values  for  the  rest  of  the  metals,  with  the  exception 
of  tungsten,  agree  fairly  well  with  those  given  in  the  literature.  For  tungsten, 
our  lower  y  j  value  undoubtedly  arises  from  our  value  of  a  which  is  lower  than 
that  ususally  found  in  the  literature.  Our  or  value  was  calculated  from  very 
recent  data.  However,  the  accuracy  of  this  data  needs  to  be  evaluated. 

B.  yQ  FROM  THERMAL  EXPANSION  DATA  AT  ZERO  PRESSURE 

a.  Calculations 

As  discussed  in  Section  4.  A.c.,  the  constancy  of  yG  with  temperature 
may  be  tested  by  the  relationship 

.  v00o 
Vg  ^00  cv 

In  this  study  the  calculation  of  yG  by  the  above  formula  was  accomplished 
by  calculating  the  ratio  Vqq/Aqq  at  room  temperature  using  the  values  of  y^, 

Cy,  and  Q  from  Tahle  5-1.  Whether  this  procedure  is  strictly  valid  needs  fur¬ 
ther  investigation.  The  values  of  O  and  Cy  at  temperatures  other  than  room 
were  calculated  in  the  same  manner  as  inTable5-l.  The  results  of  these  cal¬ 
culations  are  presented  in  Table  5-2. 


b.  Discussion  of  Results 

The  accuracy  of  the  preceding  results  suffer,  as  did  the  results  of  the 
calculation  of  y^,  from  a  lack  of  knowledge  of  precise  values  of  coefficients 
of  thermal  expansion.  However,  it  would  appear  that,  for  the  metals,  yG 
tends  to  increase  with  increase  of  temperature  while,  for  the  carbides,  yG 
tends  to  decrease. 

The  increase  in  yG  with  temperature  is  probably  not  real  because,  as 
discussed  in  Section  4.  A.c,  the  formula  yG  =  V^cr/Z^gCy  ’•  valid  only  as 
an  approximation,  especially  at  higher  temperatures. 


r>-9 


LMSC/HREC  A784868 


Table  5-2 

yG  FROM  THERMAL  EXPANSION  DATA  AT  ZERO  PRESSURE 


Material 


Af 


Be 


W 


Temp 

°K 

a  x  10+5 
(°K) " 1 

CV 

cal  . 

(gm-atm  °K) 

V00^00 
cal  (gm  .atm)-  1 

yG 

273 

7.01 

5.496 

1.73xl05 

2.20 

293 

7.01 

5.606 

2.17 

296 

7.01 

5.603 

2.17 

298 

7.01 

5.602 

2.17 

303 

7.48 

5.623 

2.31 

373 

7.48 

5.837 

2.22 

436 

8.47 

5.997 

2.44 

573 

8.47 

6.252 

2.34 

708 

9.92 

6.385 

2.69 

773 

9.92 

6.491 

2.64 

303 

3.35 

3.88 

1.46  x  105 

1.26 

373 

3.35 

4.45 

1.10 

573 

4.91 

4.97 

1.44 

673 

4.82 

5.42 

1.30 

873 

5.59 

5.78 

1.41 

1073 

5.99 

6.03 

um 

1273 

6.46 

6.28 

■S3 

298 

11.23 

5.89 

7.08  x  105 

1.35 

400 

1 1.76 

5.96 

1.40 

500 

12.28 

6.04 

1.44 

600 

12.80 

6.11 

1.48 

700 

13.32 

6.18 

1.53 

800 

13.85 

6.26 

1.57 

900 

14.37 

6.33 

1.61 

LMSC/HREC  A784868 


Table  5-2  (continued) 


Material 

Temp 

°K 

a  x  io+5 
(°K) 

CV 

cal 

( gm- atm  °K)  * 

V00^00 
cal(gm .  atm)"  ^ 

yG 

1 

14.89 

6.40 

1.65 

15.93 

6.54 

1.72 

1400 

16.98 

6.67 

1.80 

1600 

18.02 

6.81 

1.87 

1800 

19.07 

6.93 

1.95 

2000 

20.11 

7.07 

2.01 

Ta 

298 

6.018 

5.68  x 10^ 

1.70 

400 

6.198 

1.65 

700 

6.388 

1.74 

1200 

2.24 

6.516 

1.95 

1600 

2.67 

6.658 

2.28 

2000 

2. 66 

6.893 

2.19 

2400 

3.10 

7.288 

2.42 

3000 

3.66 

8.998 

2.31 

WC 

298 

1.35 

9.816 

9.82  x 105 

1.35 

400 

1.35 

11.324 

1.17 

600 

12.582 

1.17 

800 

1.35 

13.231 

1.00 

13.698 

1.08 

1200 

1.49 

14.088 

1.04 

1400 

1.94 

14.436 

1.32 

l'.OO 

1.64 

14.757 

i 

1.09  1 

1 
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Table  5-2  (concluded) 


Material  Temp  0  x  10 


(°K) 


(gm-atm  °K)~* 


voo^oo 

cal  (gm.  atm) 


6.389 

6.443 

9.321 

10.321 

10.898 

11.505 

11.836 

12.135 

12.411 
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The  tendency  for  Y  to  decrease  with  increasing  temperature  in  the 
case  of  the  ca*  hides,  however,  is  difficult  to  explain. 

It  is  reported  (Reference  18,  p.  351)  that,  for  beryllium,  y^-.,  calculated 
using  a  variation  of  the  Gruneisen  formula  used  here,  gives  an  indication  of 
becoming  negative  at  very  low  temperatures. 

C.  y2  FOR  ALUMINUM  AND  TUNGSTEN  BY  THE  METHOD  OF  BEECROFT 
AND  SWENSON  (Reference  2) 

a.  Calculations 

The  method  of  calculating  V ^  from 


was  discussed  in  Section  4.A.e.  The  method  has  been  applied  to  both 
aluminum  and  tungsten  in  this  study. 

The  validity  of  this  method  rests  on  the  temperature  independence  and 
linear  volume  dependence  of  compressibility.  For  aluminum  the  tempera¬ 
ture  and  volume  dependency  of  compressibility  was  checked  by  plotting 
compressibilities  at  different  temperatures  v£  volume  (Figure  5-2).  As 
can  he  seen  from  the  figure,  compressibility  may  be  taken  to  a  first  approxi¬ 
mation,  in  the  pressure  range  from  1  atm  to  about  30  kb  and  in  the  tempera¬ 
ture  range  20°C  -  435°C,  as  independent  of  temperature  and  a  linear 
function  of  volume. 

Tungsten  compressibility  data,  which  was  available  at  room  tempera¬ 
ture  only,  was  found  to  be  rather  unusual  in  that  it  increased  with  decreasing 
volume  (Figure  5-3).  This  behavior  cc rtainly  would  not  be  expected  for  a 
metal,  and  it  would  be  interesting  to  examine  this  behavior  further. 

The  results  of  the  calculations  for  aluminum  ate  shown  in  Figures  5-4 
through  5-8  and  for  tungsten  in  Figure  5-9. 
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Figure  5-4  -  Zero  Degree  Isotherm  for  Aluminum 
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V(P=1  atm,  T=20  C 


Variation  of  y,  as  a  Function  of  Volume  and  Temperature 
for  Aluminum4' (for  specific  temperatures  in  the  350°  -  700°K 
range,  see  Figure  5-6) 
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Figure  5-6  -  Variation  of  y ^  as  a  Function  of  Volume  and  Temperature 

for  Aluminum  -  Scale  Expanded 
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Figure  5-7  -  Thermal  Energy  as  Function  of  Temperature  and 

Volume  for  Aluminum 
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b.  Discussion  of  Results 

In  Figure  5-4  is  shown  the  calculated  zero-degree  isotherm.  It  can 
he  seen  that  the  calculated  curve  agrees  excellently  with  the  zero-degree 
values  generated  from  shock  data  (Reference  24,  p.  214). 

In  Figures  5-5  and  5-6  are  shown  the  calculated  Y ^  values  for  alumi¬ 
num  at  different  temperatures  as  a  function  of  volume.  In  Figure  5-6  the 
scale  expanded  to  show  more  detail.  The  greatest  change  in  Y ^  occurs 
below  300°K.  The  Y^  curve  at  100°K  is  very  interesting.  The  Y ^  values  at 
this  temperature  are  surprisingly  high.  In  an  attempt  to  obtain  more 
reasonable  values  at  100°K,  the  integration  interval  was  decreased  in  the 
computation  of  energy  from  the  specific  heat  curve.  However,  the  diffi¬ 
culty  persisted,  causing  us  to  examine  more  closely  the  applicability  of  .he 
basic  assumptions  of  the  method  to  metals  less  compressible  than  sodium. 

These  assumptions  are  that  the  compressibility  is  a  linear  function  of 
volume  and  independent  of  temperature.  For  aluminum,  these  conditions 
are  only  approximately  fulfilled  as  shown  in  Figure  5-2.  We  may  conclude 
then  that  the  method  of  Beecroft  and  Swenson  is  very  sensitive  to  the  slope 
of  the  compressibility  vs  volume  curve,  and  the  accuracy  of  Y?  values  calculated 
by  it  is  doubtful  in  cases  where  the  linear  volume  dependence  and  tempera¬ 
ture  independence  of  compressibility  does  not  rigorously  hold.  Thus,  the 
values  given  in  Figures  5-5,  5-6  and  5-9  for  aluminum  and  tungsten 
prohahly  are  not  very  accurate,  especially  those  for  tungsten.  For  tungsten 
compressibility  data  were  found  for  only  one  temperature,  so  that  the 
assumption  of  volume  linearity  and  temperature  independence  of  compres¬ 
sibility  is  rather  far  fetched. 

Figures  5-7  and  5-8  show  the  variation  of  the  thermal  energy  and  the 
total  pressure  with  volume  and  temperature.  They  are  included  in  this 
report  for  general  information  purposes. 
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Figure  5-9 


Y  Values  for  Tungsten  as  Calculated  by  the  Method 
2 

of  Beecroft  and  Swenson 
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D.  GRUNEiSEN'S  PARAMETER  FOR  ALUMINUM,  BERYLLIUM,  AND 
TUNGSTEN  FROM  SHOCK  DATA 


A  comparison  has  been  made  of  the  GrQneisen  parameter  as  computed 
by  several  different  methods  from  experimental  P-V  data.  These  results 
have  been  plotted  as  a  function  of  volume  ratio  in  Figures  5-10,  5-11  and 
5-12  for  aluminum,  beryllium  and  tungsten,  respectively.  All  of  the 
curves  were  plotted  from  data  presented  in  Reference  15  with  the  exception 
of  the  curve  labeled  "Walsh,  et  al"  which  was  computed  from  the  polynomial 
least  square  curve  fit  data  presented  in  Table  III  of  Reference  24. 


In  the  Beers  report  (Reference  15)  the  experimental  shock  Hugoniot 
data  (or  the  isothermal  P-V  data  in  the  case  of  Bridgeman's  data)  were  fitted 
to  two  different  mathematical  models  which  will  be  referred  to  as  the 
Murnaghan  equation  and  the  linear  equation.  The  Murnagkan  equation  has 
the  form 


(If-, 


(1) 


For  each  material,  values  of  A  and  K  were  chosen  to  minimise  the  error 
between  the  mathematical  expression  and  the  experimental  data. 


The  linear  model  assumes  a  linear  relationship  between  the  shock 
velocity  and  particle  velocity  leading  to  an  equation  of  the  form 


U 


=  c  ♦  s  u 

0  P 


(2) 


Again,  for  each  material,  values  of  and  S  were  chosen  to  minimise  the 
error  between  the  linear  model  and  the  experimental  data.  By  combining 
Equation  2  and  the  Rankine -Hugoniot  relations  the  following  equation  is 
formed.  r2  iv  -  v» 

P  .  A(W_Y).,t  ,,, 


{VOH-S,V0,fV>f 


Equations  1  and  3  represent  two  different  analytical  expressions  for 
the  same  experimental  P-V  data. 
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In  the  Beers  report,  the  Slater  formula  and  the  Dugdale- Mac  Donald 
formula  were  applied  directly  to  Equations  1  and  3  to  obtain  values  of  the  Gr&neisen 
parameter.  The  application  of  the  Slater  and  Dugdale -MacDonald  formulae 
to  Rankine-Hugoniot  and  to  non-zero  isothermal  data  is  not  theoretically 
correct,  since  the  derivations  of  both  formulas  were  based  on  isentropic 
P-V  relationships.  However,  for  very  small  compressions  the  slope  and 
curvature  of  the  Rankine-Hugoniot  curve  approach  that  of  an  isentrope,  which 
suggests  that  the  Beers  approximation  of  the  Grttneisen  parameter  may  be 
fairly  good  for  small  compressions. 

The  Gruneisen  parameter  computed  by  Walsh,  etal,  (Reference  24) 
was  obtained  by  generating  a  zero-degree  isotherm  using  experimental  shock 
data  and  commonly  available  material  properties.  The  Slater  formula  was 
then  applied  to  the  zero-degree  isotherm  to  obtain  the  Gruneisen  parameter. 
Since  the  aero-degree  isotherm  is  also  an  isentrope  this  method  is  compatible 
with  the  assumptions  inherent  in  the  Slater  formula. 


The  Los  Alamos  equation  is  based  on  the  linear  equation  and  the 
Dugdale -MacDonald  formula.  Using  Equation  3  and  the  Dugdale -MacDonald 
formula,  the  Gruneisen  parameter  and  the  slope  of  the  Gruneisen  parameter 
with  respect  to  volume  are  computed  for  zero  compression.  The  Gruneisen 
parameter  at  other  volumes  is  then  computed  by  the  following  equation 

'  ’  V<VOH-V,($p.0  (4) 


The  following  two  expressions  were  obtained  from  Reference  13. 

y0H=  2  S  "  1 


(5) 

(6) 


Substitution  of  Equations  5  and  6  into  Equation  4  leads  to  Equation  7  which 
is  called  the  Los  Alamos  equation  in  the  Beers  report. 


Y  « 


(7) 
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The  only  experimental  data  required  to  compute  the  Grunieaen  parameter 
by  thia  equation  are  the  ahock  and  particle  velocitiea  which  makea  it 
attractive  for  materiala  for  which  limited  information  ia  available. 


Equation  5  haa  been  verified  by  independent  derivation.  However,  an 
attempt  to  derive  Equation  6  reaulted  in  the  following  expreaaion: 


(8) 


Until  the  diacrepancy  between  Equation  6  and  Equation  8  ia  resolved,  the 
validity  of  Equation  7  ia  in  doubt. 


E.  r2  USING  MURNAGHAN'S  EQUATION 


One  of  the  methods  investigated  as  a  means  of  evaluating  the 
Gruneisen  parameter  involves  the  use  of  the  Murnaghan  equation  to  repre¬ 
sent  P-V-T  behavior.  The  basic  assumption  inherent  in  the  Murnaghan 
equation  is  that  the  isothermal  bulk  modulus  of  a  material  is  a  linear 
function  of  pressure  (see  Reference  26).  From  this  basic  assumption  the 
following  P-V  relationship  along  an  isothermal  can  be  derived. 


OT 

7 - 

OT 


(; 


B 


OT  \ 

r7 


OT 


-  1 


(9) 


in  which 

P  =  pressure 
V  =  volume 
B  =  bulk  modulus 

B'  =  rate  of  change  of  bulk  modulus  with  pressure 


and  the  subscript  OT  indicates  that  the  value  of  the  variable  should  be  taken 
at  0  pressure  and  the  temperature  of  the  isothermal.  Thus,  a  complete 
P-V-T  relationship  exists  for  a  substance  when  B,  B^  and  V  are  known  as 
a  function  of  temperature. 
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The  objective  of  thi6  investigation  was  to  evaluate  the  Gruneisen 
parameter,  Y^,  as  a  function  of  volume,  V,  and  temperature,  T.  The  basic 
equation  used  in  this  evaluation  is 


V(P-  P0) 

y  =  - ii. 

2  E-E0 


(10) 


in  which  E  is  energy  at  the  specified  volume  and  temperature,  and  Eq  and  Pq 
are  the  energy  and  pressure  at  zero  degrees  absolute.  At  a  given  temperature 
the  numerator  of  Equation  (10)  can  be  expressed  in  terms  of  V  and  T  by  use 
of  Equation  (9).  It  remains  to  express  (E  -  Eq)  in  terms  of  V  and  T. 


The  path  followed  (1-2- 3-4)  in  evaluating(E-EQ) and  the  subscripts  used 
in  the  following  equations  are  illustrated  in  Figure  5-13.  The  energy 
difference  between  0  and  00  is 


(ID 


From  00  to  0T  the  change  in  energy  is 


E 


0T 


C  (T)  dT 
P 


(12) 


and  from  0T  to  PT 

EpT  -  Eqt 


P(V)  dV 


But 


in  which 

or  =  coefficient  of  thermal  expansion 
B  =  isothermal  bulk  modulus 
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If  the  assumption  is  made  that  OB  is  a  function  of  temperature  only  (Refer* 
ence  2) 


a 


OT 


B 


OT 


and 


EPT  ‘  E0T 


=  T  O 


OT  B0T 


(VpT  -  V0T) 


•/ 


PT 


P(V)  dV 


(13) 


OT 


The  sum  of  Equations  (11),  (12)  and  (13)  is  the  energy  difference  E  -  Eq  in 
Equation  (10).  P(V)  along  the  two  isothermals  in  Fipire  5-13  can  be 
expressed  in  terms  of  Equation  (9)  when  the  appropriate  subscripts  are  used. 


These  equations  were  programmed  and  Y ^  was  computed  for  sodium 
and  aluminum.  Input  data  for  these  computations  are  presented  in  Table 
5-3  and  Figures  5-14  through  5-18. 


Volume  expansivity  data  are  given  in  Table  5-3  for  a  number  of 
materials.  The  numbers  appearing  in  Table  5-3  are  coefficients  deter¬ 
mined  by  least  squares  techniques  for  the  following  polynomial: 

V  -  V 

=  0.01  (A  +  BT  +  CT2  +  DT3)  (14) 

std 


in  which 

V  =  volume 

V  .  =  volume  at  293°K 
sta 

T  =  temperature,  °K 

The  data  from  which  the  coefficients  were  determined  were  obtained 

from  Reference  2  for  sodium  and  from  References  8-10  for  the  other  materials. 


In  the  computation  of  Y7  a  value  of  23.66  cm3/mole  was  used  for  the 

^  3 

standard  volume  of  sodium,  and  a  value  of  9.998  cm  /mole  was  used  for 
aluminum. 
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Figure  5-13  -  Path  Followed  in  Evaluating  Energy  Difference 
Term  in  Equation  (10) 
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Figure  5-14  -  Specific  Heat  of  Sodium  (From  Reference  14) 
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Figure  5-15  -  Specific  Heat  of  Aluminum  (From  Reference  14) 
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Figure  5-16  -  Isothermal  Bulk  Modulus  of  Sodium  (From  Reference  2) 
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Figure  5-17  -  Isothermal  Bulk  Modulus  of  Aluminum  (From  Reference  28) 


*  * 

*  t 

♦  • 

l 

1 

* 

!«« 

n 

LMSC/HREC  A784868 


Figure  5-18  -  Pressure  Derivative  of  Isothermal  Bulk  Modulus  of 
Sodium  (from  Reference  28) 
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Table  5-3 

VOLUME  EXPANSIVITY  DATA 


1 - 

Material 

A 

B  x  103 
(°K)’  1 

C  x  108 
(°Kf2 

DxlO11 

(°K)"3 

Na 

-4.14 

-1.88 

8660 

-10700 

A i 

-1.36 

2.53 

797 

297 

Be 

-1.40 

4.59 

-1.25 

92.7 

W 

-0.404 

1.28 

2.11 

4.69 

Ta 

-0.593 

2.04 

-2.17 

7.67 

WC 

-0.424 

1.43 

-13.1 

11.4 

SiC 

-0.384 

1.32 

-15.4 

-0.13 

The  pressure  derivative  of  isothermal  hulk  modulus  for  aluminum 
used  in  the  computations  has  a  constant  value  of  5.55.  This  value  was 
obtained  from  Reference  29.  Both  the  isothermal  bulk  modulus  and  the 
pressure  derivative  of  isothermal  bulk  modulus  for  aluminum  were  derived 
by  ultrasonic  techniques. 

The  computed  values  of  for  sodium  and  aluminum  are  plotted  in 
Figures  5-19  and  5-20,  respectively.  For  small  compressions  the  results 
are  similar  to  the  results  obtained  by  the  Beecroft  and  Swenson  method. 
Again,  the  curve  for  aluminum  at  100°K  is  unrealistically  high.  Of 
greater  interest  are  the  negative  values  obtained  for  at  high  compressions. 
Although  Pautamo  in  Reference  27  refers  to  negative  values  of  the  Gruneisen 
parameter,  it  seems  more  likely  in  this  case  that  the  negative  values  are 
due  to  assumptions  inherent  in  the  method  and  in  limitations  of  the  input  data. 

Pressure  curves  corresponding  to  the  Y ^  curves  are  shown  in  Figures 
5-21  and  5-22. 
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Figure  5-19  -  y,  as  a  Function  of  Volume  and  Temperature  for  Sodium 
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Figure  5-22  -  Pressure  as  a  Function  of  Volume  and  Temperature  for 
Aluminum 
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Section  6 

FOAMED  AND  DISTENDED  MATERIAL 

A.  HIGH  PRESSURE  REGION 

A  shock  adiabat  of  a  foam  suitable  for  discussion  in  the  high  pressure 
region  is  shown  by  the  dotted  line  in  Figure  6-1.  The  solid  line  represents 
the  shock  adiabat  of  solid,  compacted  material. 

The  volume  is  the  volume  of  the  foam  and  the  volume  of  non- 
porous  material.  The  volume  contains  the  same  amount  of  mass  as 
volume  V^;  the  densities,  however,  are  different.  In  other  words,  a  foam 
or  original  volume  is  crushed  to  volume  VQH  at  P  essentially  equal  to 
zero. 


Figure  6-1  -  P  v£  V  Curve  to  Illustrate  Relationship  between  Shock 
Compression  of  Compacted  Solid  and  of  Foam 
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The  different  shock  adiabats  follow  from  the  fact  that  a  shock  compres¬ 
sion  of  a  foam  is  a  more  irreversible  process  than  the  corresponding 
compression  of  a  compacted  solid. 

At  volume  V the  shocked  compacted  solid  and  the  shocked  foam  have 
the  s?.me  static  energy.  This  is  given  by  the  zero  degree  isotherm  at  volume 
.  The  difference  in  internal  energy,  therefore,  is  caused  by  a  greater 
thermal  energy  in  the  shock  foam,  i.e., 

E„  -  E/  =  AE(,  , 

H  f  thermal 

Similarly,  the  difference  in  pressure  at  volume  is  the  result  of  a 
difference  in  thermal  pressure: 

f  thermal 

The  Grdneisen  parameter  therefore  can  be  determined  from 


'  The  preceding  discussion  was  taken  mainly  from  Krupnikov,  Brazhnik, 
and  Krupuikova,  1962;  Zel'dovich  and  Rayzer,  1965,  p.  874;  and  Morgan, 
Rockowitz,  and  Atkinson,  1965. 

The  experimental  data  for  tungsten  foams  and  for  aluminum  foams 
(Morgan,  Rockowitz,  and  Atkinson,  1965,  p.  95),  can  be  represented  to  a  first 
approximation  by  the  equation 

ET  =  g(VPT)  +h(VPT)2 


6-2 


I. 

i. 

L 


LiMSC/HREC  A784868 


where  g  and  h  are  constants.  The  Gruneisen  parameter,  therefore,  can  be 
represented  as 

ldET\  1 

Y3  =  V\aPl)  =  g  +  2hVP 
1  V  1 

Because  E^,  is  a  function  of  VP^  only.y^  is  a  function  of  E^  only. 

The  Gruneisen  parameter  y ^  determined  in  this  fashion  for  aluminum 
foams  shows  a  dependence  on  the  thermal  energy.  This  is  shown  in  Figure  6-2. 

It  should  be  pointed  out  that  the  reproducibility  and  validity  of  shocked 
foam  data  is  still  questionable.  (See  Morgan,  Rockowitz  and  Atkinson,  1965, 
p.  107,  and  McQueen,  1964,  p.  69.) 


E  Thermal  Energy,  Joules/gm 


Figure  6-2  -  Gruneisen  Parameter  as  a  Function  of  Thermal  Internal 
Energy  for  Aluminum  (From  Morgan,  Rockowitz,  and 
Atkinson,  1965,  p.  104) 
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Russian  workers  (see  Al'tshuler,  1965)  have  expanded  the 
Mie-Gruneisen  equation  of  state  to  take  into  account  electronic  contributions 
at  high  temperatures  and  pressures.  Their  equation  is  essentially 


a  e 


where  the  subscripts  a  and  e  refer  to  atomic  and  electronic  contributions. 
The  Ye  parameter  apparently  varies  rather  slowly  with  density,  going  from 
l/3  for  an  isolated  atom  at  zero  pressure  to  the  limit  2/3  for  an  ideal 
Fermi-Dirac  gas  (Brush,  1967). 

In  connection  with  the  equation  of  state  in  foams,  Bjork,  et.  al.,  (1967) 

modified  existing  equations  of  state  to  agree  with  experiments  conducted  by 

the  Russians  on  porous  samples  and  to  fit  other  available  data.  In  this 

treatment  the  electronic  thermal  contribution  to  pressure  and  energy  are 

included.  The  equation  was  incorporated  into  a  hydrodynamic  code  by 

constructing  a  table  consisting  of  a  matrix  of  pressure  elements  P^,  where 

i  and  j  correspond  to  certain  values  of  rj  and  E.  The  hydrodynamic  code 

could  then  compute  the  pressures  it  needed  by  two-way  interpolation  in  the 

table.  It  appears  that  no  allowance  for  the  variation  of  Y  with  temperature 

a 

has  been  made. 

Figure  6-3  shows  a  comparison  of  PV  curves  generated  by  this 
treatment  with  experimental  data. 
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P(mbar§) 


Figure  6-3  -  Comparison  of  Aluminum  Hugoniots  with  Experimental 
Data  (Bjork,  et  al.t  1967) 
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B.  COMPLICATIONS  CAUSED  BY  STRESS  ANISOTROPY  AND  HYSTERESIS 

EFFECTS 

a.  Introduction 

The  preceding  discussion  concerned  situations  where  the  pressure 
could  be  considered  to  be  the  same  as  a  hydrostatic  pressure  and  the  states 
of  the  materials  before  and  after  shock  or  hydrostatic  compression  to  be 
those  of  thermodynamic  equilibrium.  Actually,  one  or  the  other  of  these 
two  assumptions  often  break  down  in  practice  and  may  lead  to  unpredictable 
shock  behavior  patterns.  For  example,  in  foams  and  solid  ductile  materials, 
shock  compression  at  low  pressures  may  exhibit  a  shock  front  which  splits 
into  two  fronts  (Dudley,  1966;  Rempel,  1963)  —  the  first  front  known  as  the 
elastic  precursor  and  the  second  as  a  plastic  shock  wave.  Experimental 
confirmation  of  various  mathematical  models  for  such  behavior  has  been 
good  in  some  cases  and  not  so  good  in  others. 

Attention  was  focused,  in  reviewing  these  mathematical  models,  on  how 
the  equations  of  state  were  formulated  for  foams  in  the  low  pressure  regions, 
especially  for  the  BCD  part  of  the  shock  compression  P-V  curve  which  is 
well  approximated  as  shown  in  Figure  6-4. 

It  became  apparent  almost  immediately  that  no  clear,  systematic 
treatment  of  the  basic  thermodynamics  and  physics  of  hysteresis  effects 
exists.  The  formulation  in  the  literature  of  equations  of  state  for  materials 
displaying  hysteresis  effects  are  extremely  difficult  to  understand  and  evaluate 
because  most  of  these  formulations  are  "couched  in  terms  of  elasticity  theory 
and  dislocation  concepts"  (passage  in  quotes  taken  from  Fitzgerald,  who  is 
quoted  again  in  the  next  paragraph).  The  same  criticism  may  be  made  of 
many  current  explanations  of  mechanical  behavior.  Fitzgerald  (1966)  neatly 
sums  up  the  confusion: 
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"Current  explanations  of  mechanical  properties  of  solids  are 
generally  couched  in  terms  of  elasticity  theory  and  dislocation 
concepts.  The  latter  are  employed  chiefly  in  attempts  to  ex¬ 
plain  obvious  nonelastic  behavior  such  as  plastic  deformation, 
but  even  in  this  case  elastic  ideas  prevail  to  a  surprising  ex¬ 
tent.  That  is,  calculations  of  the  forces  or  interactions  between 
dislocations  and/or  other  crystal  lattice  defects  are  often  made 
in  terms  of  linear,  elastic  continuum  surrounding  the  imperfec¬ 
tions.  Thus  a  clearly  nonelastic  phenomenon  like  plastic  slip 
is  treated  in  terms  of  linear  elastic  concepts  and  constants. 

After  thirty  years,  this  approach  has  led  to  no  numerical  calcu¬ 
lations  in  agreement  with  measured  values  of  any  nonelastic 
quantity,  except  through  .he  arbitrary  assignment  of  values  to 
one  or  more  constants,  and  even  qualitative  provisions  for  the 
existence  of  certain  well-known  classes  of  mechanical  behavior 
are  missing. " 


Fitzgerald  goes  on  to  develop  a  treatment  of  mechanical  properties  in 
terms  of  wave  mechanics,  which  leads  directly  to 

"predictions  of  nonelastic  audiofrequency  resonances  in  vibration 
experiments  and  to  anticipation  of  acoustic  emission,  simple  and 
quantitative  explanations  of  hypervelocity  impact  phenomena, 
stress -strain  relations  for  poly  and  single  crystals,  sliding 
friction  and  the  formation  of  crystal  mosaics.  Also  obtained  is 
the  correct  form  of  the  stress -strain  law  for  cubic  metals  and  its 
variation  with  single -crystal  orientation,  temperature,  parity, 
sample  size  and  various  kinds  of  sample  pretest  treatment." 


If  the  above  claims  are  really  true,  Fitzgerald's  approach  is  a  great 
advance.  Fitzgerald's  work  certainly  appears  very  promising  and  should 
be  explored  in  more  detail.  However,  a  great  need  still  exists  for  a  clari¬ 
fication  of  irreversible  processes  in  terms  of  classical  thermodynamics  and 
basic  physics.  Bridgman's  work  (1950)  in  this  direction  indicates  that  such 
an  approach  can  yield  much  insight  and  information  simply  and  directly. 


At  the  very  least,  such  an  approach  would  enable  explorations  of  the 
limits  of  applicability  of  classical  thermodynamics  and  basic  physics  and, 
thus,  provide  a  firm  basis  for  decisions  on  approximations  in  practical 
calculations  and  the  need  in  certain  cases  for  more  complicated  theories, 
such  as  Fitzgerald' 8  or  irreversible  thermodynamics.  Regardless  of  the 
light  such  an  approach  can  throw  on  irreversible  processes,  a  clarification 
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of  the  various  definitions  of  "state"  is  imperative  if  equations  of  state  are  to 
be  formulated  at  all. 

In  this  report  only  a  beginning  of  a  clarification  of  irreversible 
processes  and  definitions  of  "state"  in  terms  of  classical  thermodynamics 
and  physics  can  be  attempted  because  of  time  limitations.  However,  from 
the  following  discussion,  it  will  be  obvious  that  the  classical  approach  is 
far  from  exhausted  or  completely  understood. 

b.  Classification  of  Mechanical  Processes  and  Stresses 

First  of  all,  it  is  necessary  to  clearly  separate  elastic,  reversible 
behavior  from  irreversible,  hysteresis  type  behavior  and  internal  stresses 
from  external  stresses.  Elastic  or  reversible  displacements  are  those  in 
which  no  work  is  dissipated  as  heat.  In  an  elastic  displacement,  the  dis¬ 
placing  force  can  be  reversed  at  any  point  in  the  displacement  and  the 
material  brought  back  to  the  original  state.  All  irreversible  displace¬ 
ments  involved  work  dissipation,  and  reversal  of  the  displacement  force 
may  or  may  not  bring  the  material  back  to  its  original  state.  Also,  a  con¬ 
sideration  in  classifying  hysteresis  behavior  is  whether  the  process  approaches 
reversibility  if  done  slowly  enough.  We  believe  that  mechanical  hysteresis 
behavior  probably  can  be  classified  according  to  one  of  live  or  six  general 
classes  or  combination  thereof,  according  to  how  the  original  state  can  be 
restored  and  whether  time  is  a  factor  or  not.  However,  in  the  short  time 
that  was  available,  we  are  not  sure  that  we  have  completely  covered  every 
possibility  or  that  each  of  the  five  types  is  indeed  unique. 

The  type  of  hysteresis  most  pertinent  to  the  case  of  a  foam  collapsing, 
either  by  fracture  or  flow,  is  one  in  which  the  initial  cannot  be  restored  by 
reversal  of  displacement  force,  and  no  matter  how  slowly  the  force  is 
increased  the  collapse  can  never  be  reversible. 
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In  this  respect  the  collapse  of  a  solid  foam  has  much  in  common  with 
the  crystallization  of  a  subcooled  liquid.  A  subcooled  liquid  like  a  foam  is 
a  metastable  equilibrium  state.  If  the  subcooled  liquid  is  catalyzed  or  suddenly 
jarred,  the  solid  phase  appears  and  increases  irreversibly  until  equilibrium 
is  established.  If  a  subcooled  liquid  is  cooled  slowly,  a  critical  temperature, 
called  the  temperature  of  maximum  subcooling,  is  reached,  and  the  solid 
phase  suddenly  precipitates  and  continues  to  grow  until  equilibrium  is  estab¬ 
lished.  Although  nothing  has  been  found  in  the  literature  about  what  would 
happen  to  a  subcooled  liquid  if  the  pressure  were  gradually  increased  at 
constant  temperature,  it  is  most  probable  that  a  critical  pressure  exists, 
for  as  Bridgman  notes  (1935,  p.  9)  an  increase  of  pressure  is  equivalent  to 
a  decrease  in  temperature  in  some  respects. 

The  analogy  of  a  subcooled  liquid  and  a  solid  foam  is  further  strengthened 
by  the  fact  that  a  metal  powder  sinters  slowly  at  moderate  temperatures  but 
increases  with  a  rise  in  temperature.  At  the  Tammann  temperature  (near 
0.5  Tm)  the  sintering  rate  accelerates  markedly  (Gregg,  1961,  p.  155). 

The  thermodynamic  characterization  of  metastable  states  will  be  dis¬ 
cussed  later.  We  now  turn  our  attention  to  stress  as  a  state  variable. 

If  a  hydrostatic  pressure  P  is  applied  to  an  elastic  solid,  say  a  cube 
of  material  is  immersed  in  oil  which  is  in  turn  compressed  by  a  piston,  the 
work  done  on  the  solid  will  be  given  by  -PdV,  where  dV  is  the  change  in 
volume  of  the  solid.  The  work  will  be  -PdV,  even  if  the  cube  deforms 
anisotropically  as  long  as  the  deformation  is  elastic.  In  other  words  the 
compression  may  be  more  in  one  direction  than  in  the  others,  and  the  cube 
may  be  deformed  say  to  a  rectangular  block.  Since  the  process  is  reversible, 
the  internal  P  will  be  the  same  as  the  external  applied  P.  If  we  further 
assume  that  the  temperature  is  constant,  the  state  of  the  solid  material  may  be 
characterized  by  two  independent  state  variables  P  and  T. 
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The  change  in  internal  energy  of  the  material  as  the  result  of  the  applied 
stress 


dE  =  TdS  -  PdV1  -  rf2  dL 

Thus  the  energy  is  now  a  function  of  three  independent  state  variables. 
The  above  system  discussion  illustrates  Reiss'  clarification  of  the  differences 
between  stable  thermodynamic  equilibrium  and  metastable  equilibrium. 
According  to  Reiss,  a  system  is  in  equilibrium  if  no  discernable  changes 
occur  during  the  time  span  in  which  we  are  interested.  However,  every 
state  of  equilibrium  is  subject  to  constraints  and  the  number  of  constraints 
on  a  system  determine  the  number  of  state  variables  which  must  be  specified 
in  order  to  characterize  the  equilibrium.  Furthermore,  the  ways  in  which 
a  system  can  do  work  is  determined  by  the  number  of  constraints.  Thus,  in 
the  above  example,  we  have  three  constraints  —  the  thermostat  holding  the 
temperature  at  T,  the  piston  holding  the  pressure  at  P,  and  the  applied 
tensile  stress.  Release  of  any  one  of  these  constraints,  holding  the  other 
two  constant,  enables  the  system  to  do  a  unique  kind  of  work. 

According  to  Reiss,  the  only  difference  between  metastable  equilibrium 
and  stable  equilibrium  is  the  number  of  state  variables  involved.  This  view 
greatly  simplifies  the  understanding  of  metastable  and  stable  equilibria. 
However,  a  corresponding  clarification  of  unstable  and  neutral  equilibria 
still  remains  for  the  future. 

A  further  shortcoming  of  Reiss'  work  is  that  he  fails  to  clarify  internal 
and  external  constraints.  In  an  externally  constrained  system  at  equilibrium, 
such  as  discussed  previously,  the  internal  and  external  intensive  parameters 
are  the  same.  Thus  the  pressure  is  the  same  in  the  system  and  its  surround¬ 
ings,  as  are  temperature  and  stress. 

Internal  constraints  arise  from  activation  energy  barriers.  Edelglass 
(1966,  p.  226)  presents  a  discussion  of  activation  energies  and  their  relation¬ 
ship  to  reversible  and  irreversible  processes,  which  is  very  instructive.  A 
brief  summary  of  his  discussion  is  presented  in  the  next  paragraph. 
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Graphically,  the  free  energies  of  the  various  states  may  be  repre¬ 
sented  as  shown  in  Figure  6-6. 


* -  -■  . .  '  1  >»' 

Configuration  (x) 

Figure  6-6  -  Free  Energy _vs  Configuration  Paths  for 
Irreversible  "Reactions" 

The  states  A,  0,  and  C  are  the  metastable,  activated  and  stable 
equilibria  states.  The  activated  state  B  is  a  state  of  "unstable"  equilibria. 
The  term  "unstable"  is  put  in  quotes  because  it  is  not  clear  whether  the 
activated  unstable  state  ("The activated  complex  is  stable  for  atomic  dis¬ 
placements  in  all  directions  but  one."  —  Glasstone,  et  al. ,  p.  185)  is  the 
same  as  unstable  equilibrium  usually  referred  to  in  thermodynamic  texts. 
The  term  AF  is  the  free  energy  change  for  the  transformations  A-^-C 
or  C— ►A,  while  AF^  is  the  free  energy  of  activation  for  A— ♦■B  and 
AFq  for  C—^B.  Thus,  once  A  has  passed  over  B  and  fallen  into  state 
C,  there  is  little  chance  it  will  cross  the  barrier  B  to  go  back  to  A,  and 
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the  resultant  reaction  will  be  A  -*■  C ,  or  an  irreversible  change  at  a  rate 
which  is  determined  by  the  frequency  of  passage  over  the  barrier  B. 

Contrast  this  behavior  to  a  reversible  change,  whose  free  energy  vs. 
configuration  path  may  be  represented  as  shown  in  Figure  6-7. 


F 


■■  ^ 
Configuration  (x) 


Figure  6-7  -  Free  Energy  v£  Configuration  Paths  for 
Reversible  "Reactions" 

Here,  the  change  of  A  -►C  is  just  as  great  as  C-*-A. 

Under  stress  the  above  figure  assumes  the  form  of  the  preceding  one, 
so  that  there  is  a  net  flow  in  the  direction  of  the  force.  This  will  result  in 
a  strain  that  is  rate  controlled,  and  the  rate  will  have  an  exponential  tempera¬ 
ture  dependence.  If  the  stress  is  sufficient  to  overcome  the  barrier,  the 
strain  deformation  will  no  longer  be  temperature  dependent  (Edelglass, 
p.  228). 
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More  discussion  of  the  above  type  processes  would  point  out  many  more 
interesting  aspects  of  rate  controlled  phenomena  in  solid  behavior. 

In  connection  with  metastable  equilibria,  the  question  arises  as  to 
which  equilibria  —  stable  or  meta stable  —  will  be  more  stable  under  certain 
conditions.  This  question  arises  because  apparently  there  are  temperatures 
where  a  solid  powder,  which  we  would  classify  as  a  metastable  equilibria, 
is  more  stable  than  a  solid  crystal,  a  stable  equilibrium  case.  Huettig 
(Kingston  and  Huettig,  1951)  suggests  that  above  T  =  893°K  powdered  gold 
is  more  stable  than  solid  crystal  gold  (mp.  1337°K).  If  it  is  really  the  case, 
interesting  ramifications  may  result  in  real  situations  such  as  an  X-ray 
burst  powdering  a  missile  casing.  The  above  situation  is  perhaps  analogous 
to  the  case  of  rapid  cooling  of  liquids.  If  certain  liquids  such  as  glycerol 
are  cooled  rapidly,  they  may  be  taken  down  to  temperature  far  below  (even 
below  their  temperature  of  maximum  subcooling)  the  melting  point  of  their 
solid  phases  and  remain  in  their  metastable  states  at  that  temperature 
indefinitely  (Walton,  1966). 

The  above  brief  discussion  should  be  sufficient  to  illustrate  how  a 
unified  approach  to  hysteresis  and  rate  controlled  processes  in  terms  of 
classical  thermodynamics  and  basic  physics  is  evolving.  We  will  now  turn 
our  attention  to  explicit  thermodynamic  approaches  to  the  problem. 

c.  Unified  Thermodynamic  Approaches  to  Mechanical  Hysteresis 

In  a  very  remarkable  paper,  Bridgman  (1950)  examines  some  extensions 
of  the  methods  and  definitions  of  classical  thermodynamics  to  two  cases  of 
irreversible  mechanical  phenomena.  In  future  work  we  hope  to  go  more  fully 
into  his  reasoning  and  deductions.  Here  only  a  brief  summary  is  presented 
of  his  conclusions  about  plastic  deformation  and  processes  which  generate 
internal  disorder. 
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Consider  the  deformation  of  an  ideal  plastic  as  shown  in  Figure  6-8 
(Bridgman,  1950). 


Figure  6-8  -  The  StresB-Strain  Relations  in  the  "Ideal"  Plastic  Body- 


In  the  above  figure,  the  inclined  lines  correspond  to  elastic  deformation 
and  the  horizontal  lines  to  plastic  deformation.  For  a  body  at  point  or  "state" 
D,  the  entropy  is  the  same  as  point  A.  In  other  words,  the  entropy  in  the 
body  arising  from  plastic  deformation  is  zero. 

For  bodies  which  can  change  their  internal  order,  like  binary  alloys, 
Bridgman  introduces  a  parameter  a,  which  specifies  the  degree  of  internal 
disorder.  He  then  writes  the  entropy  of  the  system  as 

S  =  S(P,  T)  +  S(a) 

More  study  of  Bridgman's  work  should  reveal  how  the  parameter  a 
can  be  experimentally  determined  and  how  the  equation  of  state  for  a  body  at 
point  D  in  the  plastic  deformation  case  can  be  formulated  from  experimental 
stress -strain  curves. 
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In  a  group  of  papers  {Everett,  et  al. ,  1952  and  1953)  a  unified  general 
approach  to  hysteresis  is  also  approached.  However,  only  one  type  of 
hysteresis  is  dealt  with.  This  is  the  kind  where  a  stress  cycle  returns  the 
material  to  its  original  state,  but  the  paths  in  the  forward  and  reverse 
directions  cannot  be  made  to  coincide  no  matter  how  slowly  the  process  is 
carried  out.  This  type  of  hysteresis  is  a  common  phenomena  in  many  fields, 
such  as  magnetization  of  ferromagnetics,  solid  transitions  in  crystals  and 
alloys,  application  of  stress  to  solids,  adsorption  of  gases  by  solids,  etc. 
Everett  and  Smith  advance  a  domain  theory  to  explain  this  type  of  hysteresis. 
In  this  theory  the  static  macroscopic  state  of  a  system  exhibiting  hysteresis  is 
defined  at  any  instant  when  an  internal  variable  £  is  specified,  in  addition  to 
the  usual  external  independent  variables.  The  variable  £  represents  the 
fraction  of  molecules  of  the  system  in  State  II;  presumably  State  II  is  a 
meta stable  state,  the  rest  of  the  system  being  in  State  I.  Everett  and 
Smith's  treatment  will  need  to  be  studied  in  more  detail.  Their  definition 
of  metastable  equilibrium  is  confusing. 

Incidentally,  it  may  be  worth  mentioning  two  notes  of  caution  appearing 
in  the  literature  concerning  the  application  of  Onsager's  irreversible  thermo¬ 
dynamics.  Bridgman  (1950,  p.  63)  states: 

"As  far  as  1  know,  however,  all  the  examples  treated  by  the 
Belgian  school  are  such  that  some  reversible  method  exists  for 
getting  from  any  state  of  the  system  to  another,  so  that  the 
entropy  is  defined  and  may  be  evaluated  by  conventional  methods. 

The  extension  of  the  entropy  concept  to  situations  in  which  the 
entropy  may  not  be  so  defined  does  not  appear  to  have  been 
considered." 

Tisza  (1966,  p.  33)  states,  "This  note  is  merely  a  warning  that  the  very 
simple  consideration  of  irreversibility  should  not  be  thoughtlessly  applied 
to  situations  which  in  fact  involve  the  question  of  retrieval." 
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d.  Approximate  Calculation  of  an  Equation  of  State  for  a  Partially 

Collapsed  Foam 

A  calculation  of  a  Y  parameter  for  foams  based  on  a  rigorous  applica¬ 
tion  of  thermodynamics  would  obviously  involve  much  more  time  than  is 
available  in  the  present  study.  The  preceding  discussions,  however,  do  suggest 
that  the  problem  can  perhaps  lend  itself  to  some  approximations.  For  example, 
consider  the  collapsing  part  of  the  P-V  curve  for  a  foam 


where  V  and  represent  the  volumes  of  a  given  mass  of  compacted  solid 
and  of  foamed  solid  (P  =  0,  T  =  initial).  State  3  represents  the  equilibrium 
end-point  state  after  application  of  a  pressure  Pc»  If  P£  is  released 
slowly  and  adiabatically  state  3  would  relax  isentropically  to  state  4  (assump¬ 
tion).  Let  us  assume  that  state  3  represents  a  state  containing  (Vg-Vi/fVf-V) 
parts  foamed  solid  and  the  rest  compacted  solid.  Let  us  further  assume  that 
state  3,  on  cooling  to  room  temperature  at  constant  volume,  follows  the 
reversible  path  3-5.  A  7  parameter  may  then  be  given  by 
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so  that 


P 

c 


o 


or  assuming  7  is  independent  of  pressure 


c 


If  it  is  further  assumed,  as  the  Russians  did  (Krupnikov,  et  al. ,  1962),  that 
the  energy  of  a  foam  at  initial  conditions  is  little  different  from  that  of  the 
compacted  solid,  then  E^  can  be  taken  as  the  energy  of  compacted  solid  at 
the  initial  conditions. 


If  E^  can  be  determined,  and  it  probably  can  be  approximated  if  heat 
capacity  data  on  the  compacted  solid  and  foam  are  available,  then  7  can  be 
calculated. 
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Section  7 

CONCLUSIONS  AND  RECOMMENDATIONS 

Examination  of  the  theoretical  bases  of  the  Mie -Gruneisen  equation  of 
state  and  the  results  of  calculations  of  the  various  Gruneisen  parameters 
have  shown  that  its  validity  cannot  be  assumed  lightly  under  all  conditions. 

The  study  has  also  shown  that  a  knowledge  of  the  temperature  and 
volume  dependence  of  y  ^  would  not  only  provide  a  firm  basis  on  which  to 
judge  the  validity  of  the  Mie -Gruneisen  state  equation  but  also  would  permit 
direct  use  of  the  Mie -Gruneisen  type  equation  of  state 

y2(v.T) 

P-P0  =  —  y—  (E-E0) 


in  shock  wave  calculations. 

Thus,  it  would  be  most  desirable  that  methods  of  obtaining  y^  from 
experimental  data,  both  under  hydrostatic  and  shock  conditions,  be  fully 
explored.  In  particular,  it  is  suggested  that  the  accuracy  of  the  method 
of  calculating  y  ^  originated  in  this  study  be  evaluated  and  that  y ^  param¬ 
eters  be  calculated  for  a  variety  of  materials. 

It  is  also  recommended  that  the  physical  and  thermodynamic  bases  of 
formulating  equations  of  state  be  examined  more  fully.  Such  an  examination 
should  lead  to  a  theoretically  sound  approach  to  formulating  equations  of 
state  for  materials  characterized  by  more  than  two  state  variables,  such  as 
foams  and  distended  materials. 

Other  desirable,  work  indicated  by  the  study  is  (a)  obtaining  more  pre¬ 
cise  values  of  coefficients  of  thermal  expansion  from  percent  thermal  data 
in  the  literature,  (b)  a  concentrated  examination  of  the  factors  involved  in 
obtaining  reliable  y  values  for  plastics,  (c)  a  fuller  understanding  of  directional 
Gruneisen  parameters,  (d)  a  fuller  understanding  of  the  role  of  the  Gruneisen 
parameter  in  phase-change  (melting)  phenomena,  and  (3)  a  fuller  consideration, 
of  the  role  of  electrons  in  the  thermodynamics  of  solids. 
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